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1  Summary 

Image  processing  typically  involves  the  processing  of  large  quantities  of  data  by  computa¬ 
tionally  intensive  and  mathematically  complex  numerical  procedures.  A  key  task  in  image 
processing  is  to  search  for  transforms  that  faithfully  and  efficiently  represent  image  data. 
The  choice  of  a  transform  in  a  specific  application  can  depend  on  a  wide  range  of  possi¬ 
bly  conflicting  criteria,  including  the  nature  of  the  redundancies  and  noise  introduced  by 
sensors  and  image  propagation  medium  characteristics,  the  overall  accuracy,  stability  and 
computational  cost  of  utilizing  the  transform  and  the  underlying  processing  strategy. 

This  report  addresses  the  need  to  design  and  implemented  application  specific  algorithms 
for  automated  target  detection  and  classification  by  providing  novel  image  transforms  from 
nonabelian  group  harmonic  analysis.  Our  methods  include  procedures  for  incorporating 
prior  knowledge  of  target  size  and  shape  as  inputs  to  thresh-holding  and  filtering  of  spectral 
coefficients  defined  by  unitary  operators  associated  to  nonabelian  groups.  By  operating  on 
spectral  coefficients,  these  operators  can  be  targeted  to  specific  tasks  at  a  level  of  detail 
not  available  in  traditional  transform  methods  including  Fourier  methods,  discrete  cosine 
transform  methods  and  more  recently  developed  wavelet  methods.  When  combined  with 
a  corresponding  analysis  of  target  signal’s  unique  structures  and  textural  properties,  our 
methods  provide  powerful  tools  for  detection  and  classification  and  protection  from  false 
detection. 

We  have  investigated  the  use  of  group  theory  for  producing  transformations  that  will 
incorporate  several  classes  of  prior  image  information.  Identification  of  image  properties 
which  can  be  addressed  in  a  group  theoretic  framework  is  a  part  of  our  ongoing  research. 

Our  efforts  have  been  focused  on  constructing  transforms  for  discriminating  particular 
geometric  properties  in  images  as  well  as  locating  the  occurrences  of  such  geometry.  Encoding 
of  prescribed  pattern  into  an  invariant  subspace  is  a  design  process,  based  on  selection 
of  a  nonabelian  group.  Once  this  encoding  is  completed  however,  locating  the  pattern  is 
automatic.  Compared  to  methods  based  on  matched  filter  processing,  the  nonabelian  group 
projection  methods  will  be  global,  and  can  encode  multi-dimensional  geometric  information 
that  are  independent  of  precise  pixel  values,  or  pixel-to-pixel  variances. 

For  applications  to  image  processing,  especially  for  composite,  iterative  processing,  the 
software  library  of  the  unitary  transforms  of  complete  systems  of  orthogonal  idempotents 
and  associated  projection  operators  are  being  extended  and  structured  based  on  both  the 
image  processing  characteristics  and  group  theoretic  properties  of  the  projection  operators. 
These  operators  are  tested  extensively  for  performance  in  severe  additive  and  multiplicative 


noises. 
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2  Introduction 

During  the  last  ten  years  considerable  effort  has  taken  place  to  extend  the  range  of  ap¬ 
plicability  of  nonabelian  group  methods  to  the  design  of  new  filters  and  spectral  analysis 
methodologies  [8,  9,  5],  as  an  image  processing  tool  [11,  13,  19,  20]. 

Many  efforts  at  finding  a  significant  role  for  nonabelian  group  theory  in  DSP  and  image 
processing  applications  as  well  as  in  coding  and  communication  theory  have  been  limited  by 
some  or  all  of  the  following. 

•  The  choice  of  an  appropriate  group  or  groups  in  a  given  application  is  not  obvious. 

•  The  lack  of  a  large  class  of  groups  whose  harmonic  analysis  is  sufficiently  understood 
for  meaningful  applications. 

•  The  absence  of  a  conceptual  framework  which  relates  group  harmonic  analysis  to  phys¬ 
ically  interpretable  results. 

•  Fast  algorithms  may  exist  but  are  difficult  to  code  as  they  do  not  easily  lend  themselves 
to  modification. 

•  The  lack  of  relevant  models  for  applications  which  are  not  only  a  rephrasing  of  known 
methods  or  what  is  immediately  available. 

The  language  of  group  algebras  over  the  complex  field  will  be  used  throughout.  Some 
readers  may  be  more  familiar  with  the  identification  of  the  group  algebra  CG  with  the 
space  of  all  complex  valued  functions  on  G  under  G-convolution.  Harmonic  analysis  over  G 
usually  includes  a  description  of  the  (left)  G-invariant  and  irreducible  G-invariant  subspaces 
and  direct  sum  decompositions  of  CG  into  irreducible  G-invariant  subspaces. 

For  an  abelian  group  A,  complete  answers  to  these  problems  can  be  given  in  terms  of 
the  characters  of  A.  Moreover,  these  descriptions  are  especially  simple.  For  example  the 
irreducible  A-invariant  subspaces  coincide  with  the  one-dimensional  subspaces  spanned  by 
the  characters.  The  characters  of  A  determine  an  orthogonal  basis  of  C A  and  the  finite 
Fourier  transform  relates  the  delta  basis  A  of  C A  with  the  basis  of  characters  of  C  A.  This 
is  an  extremely  nice  answer  for  applications  since  the  characters  of  A  can  be  physically 
interpreted  as  frequencies. 

For  finite  groups  G  —  A  ^  B,  the  semidirect  product  of  normal  abelian  subgroup  A 
with  an  abelian  subgroup  B ,  the  characters  of  A  and  B  play  an  equally  important  but  more 
complicated  role  in  the  harmonic  analysis  of  A  <$  B.  We  have  developed  harmonic  analysis  of 
such  groups  and  directly  linked  this  harmonic  analysis  to  familiar  DSP  and  image  processing 
concepts.  Moreover  we  have  developed  for  groups  G  =  A  B 

•  algorithms  for  computing  bases  of  irreducible  G-invariant  subspaces  and  bases  compat¬ 
ible  with  direct  sum  decompositions  of  CG  into  G-invariant  and  irreducible  G-invariant 
subspaces  (spectral  bases). 

•  fast  algorithms  for  representing  data,  i.e.,  elements  in  CG  over  spectral  bases. 
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•  an  extensive  new  class  of  fast  unitary  transforms  for  data  analysis. 

The  direct  sum  decompositions  of  C G  studied  in  this  work  are  closely  related  to  but  are 
significantly  different  from  those  leading  up  to  the  Fourier  transform  of  C G  [10,  26].  Our 
decompositions  are  finer  but  are  not  generally  uniquely  determined. 

We  have  extended  our  computational  nonabelian  group  harmonic  analysis  technology 
to  include  semidirect  products  of  the  form  G  =  A$H,  where  A  is  an  abelian  group  and 
H  =  B  <£  C  is  the  semidirect  product  of  abelian  groups  B  and  C.  In  previous  efforts, 
complete  systems  of  primitive,  pairwise  orthogonal  idempotents  were  built  from  characters 
of  abelian  groups.  The  current  effort  is  based  on 

•  adapting  the  Mackey- Wigner  little  group  theorem  to  an  algorithm  for  constructing  such 
systems  for  general  semidirect  products  A  H,  where  A  is  abelian  and  H  is  arbitrary. 
This  algorithm  depends  on  knowing  such  systems  for  certain  subgroups  of  H. 

•  extending  the  results  on  abelian  group  characters  to  general  one-dimensional  characters 
of  nonabelian  groups. 

The  complete  algorithm  computes  systems  of  primitive,  pairwise  orthogonal  idempotents 
for  the  groups  G  under  study  in  terms  of  extensions  of  abelian  group  characters  to  one¬ 
dimensional  characters  of  subgroups  of  G. 

For  a  finite  group  G  of  the  form  A  (B  $  C),  where  A,  B  and  C  are  abelian  groups,  we 
have  developed  algorithms  for  constructing 

•  a  complete  primitive  idempotent  system  in  KG. 

•  a  basis  for  each  of  the  irreducible  G-invariant  subspaces  generated  by  the  primitive 
idempotents  of  the  system. 

•  the  expansion  coefficients  of  a  e  KG  relative  to  the  resulting  basis  of  KG. 

Algorithms  for  convolutions  are  investigated  based  on  Wedderburn  structure  theorem. 
Compared  to  previous  efforts,  convolutions  will  be  described  in  terms  of  two-sided  G- 
invariant  subspaces  of  C G. 

We  have  generalized  the  computational  nonabelian  group  harmonic  analysis  technology 
by  developing  the  theory  and  algorithms  for  abelian  extensions.  The  previous  efforts  have 
been  based  on 

•  adapting  the  Mackey-Wigner  (M-W)  little  group  theorem  to  an  algorithm  for  con¬ 
structing  such  systems  for  general  semidirect  products  A  where  A  is  abelian  and 
H  is  arbitrary.  This  algorithm  depends  on  knowing  such  systems  for  certain  subgroups 
of  H. 

•  extending  the  results  on  abelian  group  characters  to  general  one-dimensional  characters 
of  nonabelian  groups. 
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The  complete  algorithm  computes  systems  of  primitive,  pairwise  orthogonal  idempotents 
for  the  groups  G  under  study  in  terms  of  extensions  of  abelian  group  .characters  to  one¬ 
dimensional  characters  of  subgroups  of  G.  A  theorem  of  Wigner  and  Mackey  provides  a 
framework  for  constructing  complete  systems  of  primitive,  pairwise  orthogonal  idempotents 
for  semidirect  product  groups.  However,  M-W  little  group  theorem  requires  knowledge 
of  complete  systems  of  primitive,  pairwise  orthogonal  idempotents  for  certain  subgroups 
of  H  even  for  semidirect  product  groups.  The  framework  of  abelian  extensions  simplifies 
algorithms,  as  well  as  includes  groups  not  covered  by  the  previous  methods. 
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3  Mathematical  Results 

3.1  Introduction 

Prom  a  mathematical  perspective,  a  significant  part  of  classical  digital  signal  processing 
(DSP)  can  be  viewed  as  topics  in  finite  abelian  group  harmonic  analysis  [1,  6,  16,  4].  Fun¬ 
damental  DSP  operations  such  as  convolution  and  the  Fourier  transform  can  be  identified 
with  group  algebra  multiplication  and  group  algebra  direct  sum  decompositions  into  irre¬ 
ducible  group  invariant  subspaces.  This  interplay,  often  implicit,  has  been  responsible  for 
fast  algorithms  such  as  the  FFT  [7],  the  use  of  FFT  in  computing  large  size  convolutions 
and  correlations,  and  more  recently  for  the  development  of  polynomial  transforms  [27,  25] 
for  computing  convolutions. 

However  pleasing,  at  least  to  the  mathematically  inclined,  and  useful  this  group  theoretic 
approach  to  DSP,  it  is  based  on  a  seemingly  magical  relationship  between  DSP  applications 
and  finite  abelian  group  harmonic  analysis  and  the  simplicity  of  this  harmonic  analysis.  This 
magic  can  be  explained  by  the  powerful  role  played  by  abelian  group  characters  in  providing 
all  that  is  required  for  finite  abelian  group  harmonic  analysis  and  the  physical  interpretation 
of  these  characters  as  frequency. 

For  some  mathematicians,  this  group  interpretation  of  DSP  has  raised  the  potential  of 
an  equally  important  application  of  nonabelian  group  harmonic  analysis  to  DSP,  especially 
in  the  construction  of  group  transforms  and  group  filters  generalizing  the  classical  Fourier 
transform  and  convolutional  filters.  The  works  of  R.  Holmes  [14,  15],  M.  Karpovski  and  E. 
Trachtenberg  [18]  are  the  basis  of  much  of  the  research  in  this  direction.  Similar  ideas  in 
coding  theory  have  been  introduced  by  F.J.  MacWilliams  [21].  These  efforts  have  shown  some 
promise  but  for  the  most  part  are  more  interesting  to  mathematicians  than  DSP  engineers. 
A  more  successful  application  has  been  to  fast  algorithm  design  [2,  3,  10,  22,  23]. 

During  the  last  ten  years  considerable  effort  has  taken  place  to  extend  the  success  and 
range  of  applicability  of  nonabelian  group  methods  to  the  design  of  new  filters  and  spectral 
analysis  methodologies  [5,  8,  9],  as  an  image  processing  tool  [11,  13, 19,  20]  and  most  recently 
as  an  image  processing  tool  combined  with  graph  theoretic  modeling  of  image  data  [12,  24]. 

Many  efforts  at  finding  a  significant  role  for  nonabelian  group  theory  in  DSP  and  imaging 
applications  as  well  as  in  coding  and  communication  theory  have  been  limited  by  some  or 
all  of  the  following. 

•  The  choice  of  an  appropriate  group  or  groups  in  a  given  application  is  not  obvious. 
Often  the  groups  considered  are  those  which  are  best  known  to  the  researcher. 

•  The  lack  of  a  large  class  of  groups  whose  harmonic  analysis  is  sufficiently  understood 
for  meaningful  applications.  Often  the  dihedral  group  with  or  without  justification  is 
the  test  example. 

•  The  need  to  develop  a  conceptual  framework  which  relates  group  harmonic  analysis  to 
physically  interpretable  results. 

•  Fast  algorithms  may  exist  but  are  difficult  to  code  as  they  do  not  easily  lend  themselves 
to  modification. 
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•  Relevant  models  for  applications  which  are  not  only  a  rephrasing  of  known  methods 
or  what  is  immediately  available. 

The  finite  nonabelian  groups  in  this  study  have  the  form  G  =  A  <$B}  the  semidirect 
product  of  a  normal  abelian  group  A  with  an  abelian  group  B.  We  will  extend  the  results 
to  finite  groups  A  $  H,  where  H  is  an  arbitrary  group  whose  harmonic  analysis  is  known. 
This  includes  finite  groups  of  the  form  A  (B  C) ,  where  A,  B  and  C  are  abelian  groups 

and  in  particular  all  crystallographic  groups. 

In  the  applications  part,  we  will  see  how  these  nonabelian  groups  provide  for  new  classes 

of  imaging  and  DSP  models. 

The  language  of  group  algebras  usually  over  the  complex  field  will  be  used  throughout. 
Some  readers  may  be  more  familiar  with  the  identification  of  the  group  algebra  CG  with  the 
space  of  all  complex  valued  functions  on  G  under  (7-convolution.  Harmonic  analysis  over  G 
usually  includes  a  description  of  the  (left)  (7-invariant  and  irreducible  (7-invariant  subspaces 
and  direct  sum  decompositions  of  C(7  into  irreducible  G- invariant  subspaces. 

For  an  abelian  group  A,  complete  answers  to  these  problems  can  be  given  in  terms  of 
the  characters  of  A.  Moreover,  these  descriptions  are  especially  simple.  For  example  the 
irreducible  A-invariant  subspaces  coincide  with  the  one-dimensional  subspaces  spanned  by 
the  characters.  The  characters  of  A  determine  an  orthogonal  basis  of  C A  and  the  finite 
Fourier  transform  relates  the  delta  basis  A  of  GA  with  the  basis  of  characters  of  A.  This 
is  an  extremely  nice  answer  for  applications  since  the  characters  of  A  can  be  physically 
interpreted  as  frequencies. 

For  finite  groups  A  B,  the  characters  of  A  and  B  play  an  equally  important  but  more 
complicated  role  in  the  harmonic  analysis  of  A  B.  In  this  section,  we  will  present  the 
mathematical  basis  for  G  =  A  $  B  and  extensions  to  A  4  (B  $  C)  and  beyond  for  developing 

•  algorithms  for  computing  bases  of  irreducible  G-invariant  subspaces  and  bases  compat¬ 
ible  with  direct  sum  decompositions  of  CG  into  G-invariant  and  irreducible  G-invariant 
subspaces  (spectral  bases). 

•  fast  algorithms  for  representing  data,  i.e.,  elements  in  CG  over  spectral  bases. 

•  an  extensive  new  class  of  fast  unitary  transforms  for  data  analysis. 

We  will  have  as  a  result  a  large  class  of  nonabelian  groups  with  the  potential  of  wide 
applicability  to  DSP  and  image  processing  whose  harmonic  analysis  is  known  in  detail.  This 
harmonic  analysis  is  based  on  abelian  group  characters  which  directly  link  this  harmonic 
analysis  to  familiar  DSP  and  imaging  concepts. 

The  direct  sum  decompositions  of  CG  studied  in  this  work  are  closely  related  to  but  are 
significantly  different  from  those  leading  up  to  the  Fourier  transform  of  CG  [10,  26].  Our 
decompositions  are  finer  but  are  not  generally  uniquely  determined. 

Throughout  this  section,  A  and  B  denote  finite  abelian  groups  with  composition  given 
by  multiplication.  The  identity  is  always  denoted  by  1  and  the  inverse  of  x  £  A  by  x  .  The 
order  of  a  set  is  the  number  of  elements  in  the  set. 
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For  much  of  the  general  theory  K  denotes  a  field,  finite  or  infinite,  whose  characteristic 
does  not  divide  the  orders  of  A  and  B.  However  for  DSP  and  imaging  applications,  K  is  the 
field  of  complex  numbers. 

3.2  Abelian  group  A 

3.2.1  Characters  of  Finite  Abelian  Groups 

Suppose  A  is  a  finite  abelian  group  of  order  N.  A  mapping  r  :  A  — >  C*  is  called  a  character 
of  A  if  r  is  a  homomorphism  of  the  group  A  into  the  multiplicative  group  C*  of  nonzero 
elements  in  C, 

r(xy)  =  r(x)r(y),  x,  y  €  A. 

Denote  by  A*  the  set  of  all  characters  of  A. 

Denote  by  UN( C)  the  multiplicative  group  of  all  N- th  roots  of  unity  in  C.  UN(C)  is  a 
cyclic  subgroup  of  Cx  of  order  N.  Since  every  x  G  A  satisfies  xN  =  1,  r( x)  G  Un(C),  for 
every  character  r  of  A. 

In  the  following  examples  Cjv  is  the  cyclic  group  of  order  N  having  generator  x. 
Example  1  There  exists  N  characters  of  CN  over  C  defined  by 

rn(x)  =  e27ri^ ,  0  <n<N. 

Consider  the  direct  product  CNl  x  CNi  of  cyclic  groups  CNl  and  CNi.  Denoting  generators 
of  CNl  and  CN2  by  xx  and  x2,  every  element  in  CNl  x  CV2  can  be  written  uniquely  in  the 
form 

x^x^2,  0  <  nx  <  JVi,  0  <  n2  <  N2. 

The  characters  of  C^x  x  C^2  are  completely  determined  by  the  characters  of  Cpjx  and  C^2  ■ 
If  r  is  a  character  of  Cnx  x  Cn2,  we  can  define  characters  T\  and  r2  of  CW x  and  Cn2  by 

Ti(xi)  =  r(x),  r2(  x2)  =  t(x2). 

t  is  completely  determined  by  T\  and  r2  by 

r(x”lx ”2)  =  Ti(aq)niT2(a:2)n2,  0  <  n\  <  Ni,  0  <n2  <  N2, 

and  we  can  write  t  =  T\  ®  r2.  Conversely,  if  T\  and  r2  are  characters  of  Cnx  and  C^2,  then 
n  ®  r2  is  a  character  of  Cnx  x  Cjv2. 

Example  2  The  NiN2  characters 

rn,  n  =  (ni,n2),  0  <  ni  <  Nit  0  <n2<  N2, 

of  C/v,  x  Cpj2  are  defined  by 

n_-Tllml  n~’n2rn2 

Tnix^'x™2)  =  e  e  "2  ,  0  <  mi  <  Nu  0  <  m2  <  N2. 
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The  results  of  this  section  easily  extend  to  an  arbitrary  finite  number  of  cyclic  groups 


A  =  CNl  x  •  •  •  x  CNr. 


Each  character  r  of  A  is  uniquely  represented  by 


r  =  Ti  <g>  •  •  •  ®  tr, 

where  rr  is  a  character  of  CNr,l<r<  R.  Since  every  finite  abelian  group  A  is  the  product 
of  a  finite  number  of  finite  cyclic  groups,  this  completes  the  description  of  A*  for  every  finite 
abelian  group  A. 

3.2.2  Group  algebra  of  A 

The  group  algebra  CA  of  A  over  C  is  the  C-vector  space  of  all  formal  sums 

/  =  X  f(x)x *  /(x)  e  c> 

x£A 

with  C-algebra  multiplication 

f9  =  X  (X  /(*)0(*-1v))  V’  ft  9  ^  CA- 

yGA  Va:€j4  / 

Denote  by  L(A\  C),  the  C-vector  space  of  all  C-valued  functions  on  A.  Every  /  G  L(A  : 
C)  defines  a  formal  sum  in  CA 

f  =  X  /(*)* 

x£A 

and  we  can  identify  the  C-vector  space  CA  with  the  C-vector  space  L(A\  C).  The  multipli¬ 
cation  /  •  g  in  CA  corresponds  to  the  standard  convolution  in  L{A\  C) 

/  *  g{y )  =  X)  y  e  A>  ft  9  ^  L(A't  C). 

x£A 

Under  this  identification  the  delta  function  8y,  y  E  A  corresponds  to  the  formal  sum  in  CA 

X  6v(x)x 

xGA 

which  we  denote  by  y.  In  this  way,  we  can  view  A  as  a  subset  of  CA.  A  is  a  basis  of  the 
C-vector  space  CA  and  CA  has  dimension  N,  the  order  of  A.  The  identity  1  in  A  is  the 
identity  of  the  C-algebra  CA. 

We  will  usually  use  the  term  basis  for  any  subset  of  vector  space  which  if  ordered  is  a 
basis  in  the  usual  sense.  In  summation  expressions,  there  is  no  loss  in  doing  so.  However 
whenever  matrices  are  involved,  an  ordering  must  be  specified  even  if  implicitly. 

For  y  E  A  and  /  E  CA 

yf  =  X  f(x)yx  =  X  f(y~lx)x ■ 

x€A  x€A 
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We  call  yf  the  translation  of  /  by  y.  The  term  f(y~lx)  in  the  right-hand  summation 
identifies  yf  with  the  usual  definition  of  translation  in  engineering  terminology.  For  y  e  A, 
the  operator  L(y)  of  CA  defined  by 

L(y)f  =  yf,  /GCA, 

is  a  linear  isomorphism  of  the  C-vector  space  CA. 

A  subspace  V  of  the  C-vector  space  CA  is  called  A-invariant  if  for  each  y  £  A 

yV={yf:feV}c  V. 

For  g  e  CA  define  the  operator  L(g)  of  CA  by 

L(9)f  =  9f »  /  e  CA. 

Since  .  . 

L{g)f  =  gf  =  Mivf)  =  12  9(y)L(y)f  =  (l2  MUv) )  / 

yeA  \y€-4  / 

we  have  that 

Us)  =  MUv) 

y€A 

linearly  extends  the  domain  of  definition  of  L  from  A  to  CA.  In  particular,  the  two  definitions 
coincide  on  A. 

L(g),  g  E  CA,  is  a  homomorphism  of  the  C-vector  space  CA  but  is  not  necessarily  an 
isomorphism  as  is  the  case  with  L(y),  y  G  A.  We  can  have  gf  =  0  for  / ,  g  E  CA  with  both 
/  and  g  not  zero. 

For  an  A-invariant  subspace  V  of  CA  and  /  G  V,  we  have  for  all  g  €  CA, 

ms = e  £  y 

and  V  is  CA-invariant.  The  equivalence  of  A-invariance  and  CA-invariance  will  be  used 
throughout  this  work. 

3.2.3  Fourier  transform  over  A 

Decompositions  of  CA  into  direct  sums  of  A-invariant  subspaces  play  two  important  roles: 
algorithm  design  for  computing  products  in  CA  and  spectral  analysis.  The  Fourier  transform 
over  A  is  the  simplest  and  most  frequently  occurring  example. 

A  character  r  of  A  determines  the  formal  sum  in  CA,  r  =  £i€/1  t(x)x.  Multiplication 
of  characters  will  always  be  taken  in  CA  with  the  warning  that  in  many  places,  another 
multiplication  is  defined  under  which  A*  is  a  group.  The  importance  of  characters  in  studying 
A-invariant  subspaces  of  CA  is  contained  in  the  following  result. 
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Theorem  1  For  f  G  C  A  and  t  a  character  of  A, 

fr  =  /(r)r, 

where  /(r)  G  C  is  given  by 

hT)  =  E/^My-1)- 

y€A 

Proof  For  y  G  A 

yr  =  £  r(x)yx  =  E  r(y_1x)a;  =  r(y_1)  E  t(x)s  =  r(y-l)r. 


x£A  x£A 

The  theorem  follows  from 


x£A 


/r  =  E  f(y)(vT)  =  ( E  f(y)T(y  x) ) r- 

yeA  \yeA  ) 


For  two  characters  r  and  A  of  A,  since  rA  =  At,  we  have 

r  A  =  aA  =  ar  =  At, 

where  a  =  r(A)  =  A (r).  If  t  ^  A,  then  a;  =  0  and  tA  =  0.  If  t  =  A,  then 

a  =  E  T0/)r0/~1)  =  E 1  =  N 

and  r2  =  TVr,  proving  the  following. 

Corollary  1  For  two  characters  r  and  A  o/  A, 

Nr,  r  =  A, 


rA  = 


0,  r/  A. 


Corollary  2  A*  is  a  linearly  independent  subset  in  the  C-vector  space  CA. 
Proof  Suppose 

o  =  E  a(r)r>  Q:(r) e  c. 

r£A* 


By  corollary  1,  for  any  A  G  A* 

A  E  aiT)T  —  Na(\)\  =  0 

t€A* 

and  a(A)  =  0.  Since  this  holds  for  any  A  G  A*,  the  corollary  follows. 


By  Theorem  1,  the  C-subspace  spanned  by  a  character  r  is  ^-invariant, 

CAr  =  C  r. 

By  corollary  2,  we  have  the  direct  sum  of  one-dimensional  A-invariant  subspaces  of  CA, 

E©Ct. 
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3.2.4  Convolution  theorem 

By  corollary  2,  A*  is  a  basis  of  the  C-vector  space  C A  and  we  have  the  following  result. 
Theorem  2  C A  is  the  direct  sum  of  one- dimensional  A-invariant  subspaces 

ca=Y  ®Cr- 

rCAr 


By  theorem  2,  __ 

1  =  Y  a(T)T>  a(T)  e  c. 

tGA" 

For  any  X  E  A*, 

A  =  A  •  1  =  Y  a(T)^T  =  Na(X)\ 

t€A* 

and  a(A)  =  jj,  proving  the  following-. 


Corollary  3 


*!>• 

t€.A* 


By  corollary  3, 

/  =  =  /(r)r 
JV  t£A*  1  tGA* 

is  the  expansion  of  /  over  the  basis  A*  of  CA .  The  coefficient  set  of  this  expansion  (up  to 
scale  multiple  jj) 

f{r ),  T  E  A*, 

is  called  the  Fourier  transform  of  /  in  CA.  By  Theorem  1, 

f{r)  =  £/(y)r(y-1)- 

y€A 

The  basis  A*  of  CA  is  a  diagonalizing  basis  for  the  operators  L(g ),  g  E  C A,  since 

L(g)r  =  gr  —  ^(r)r,  t  E  A* 
and  a  diagonalizing  basis  for  multiplication  in  CA. 

Theorem  3  For  f  and  g  in  Cy4, 

fg  =  TfH  f(T)g{T)T. 

t€A+ 
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"  - 

=  pEE  />»(a)ta 

iV  r€i4‘  A€A* 

=  Tf  E  f(T)9(T)P 

JV  reA* 

completing  the  proof. 

An  A-invariant  subspace  W  of  C  A  is  called  irreducible  if  the  only  A-invariant  subspaces 
of  W  are  (0)  and  W.  For  r  £  A*,  Cr  is  an  irreducible  A-invariant  subspace  of  C A. 

Theorem  4  IfWis  an  A-invariant  subspace  of  C A,  then 

W=Y1  ®Ct  =  C>le> 

t€A 

where  A  =  VF  fl  A*  and  e  =  jj  SreA  r • 

Proof  For  /  £  W  and  A  £  A* 

A/  =  ^  £  Hr) Ar  =  ^/(A)A. 

JV  t€A* 

Since  VF  is  A-invariant, 

/(A) A  £  VF,  A  £  A*. 

jf  /(A)  ^  0,  then  A  £  VF  proving  that  every  /  £  VF  is  contained  in  the  C-linear  span  of  the 

set  of  characters  contained  in  VF,  proving  FF  =  £t€a  ©Ct • 

Since  e  £  VF  and  VF  is  A-invariant,  CAe  C  W.  For  any  A  £  A,  Ae  =  ^A  £  CAe,  proving 
VF  C  CAe,  completing  the  proof  of  the  theorem. 

The  factor  jj  in  the  definition  of  e  has  been  chosen  so  that  e2  =  e,  the  relevance  of  which 
will  be  made  clear  in  the  next  section. 

If  VF  is  an  irreducible  A-invariant  subspace  of  CA,  then  by  theorem  4,  VF  contains  a 
unique  r  £  A*  and  VF  =  Cr,  proving  the  following  corollary. 

Corollary  4  The  irreducible  A-invariant  subspaces  of  CA  are  given  by  Cr,  r  £  A  . 

Denote  the  complement  of  A  =  VF  n  A*  in  A*  by  Ac.  By  theorem  2  and  theorem  4,  we 
have  the  following  result. 

Corollary  5  CA  is  the  direct  sum  of  A-invariant  subspaces, 

CA  =  VF  ©  VF', 

where 

W'  =  Y,  ©Cr  =  CAe'> 

r€  Ac 


and  e  ^  X^re Ac 
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3.3  Finite  nonabelian  group  G 

3.3.1  Group  Algebras  of  Finite  Group  G 

Suppose  G  is  an  arbitrary  finite  group  of  order  N.  The  group  algebra  C  G  of  G  over  C  is  the 
C- vector  space  of  all  formal  sums 

/  =  /(*)*>  /(*)  e  c, 

tG:G 


with  C-algebra  multiplication 

$9  =  E  fe  /(u)2(u-1*))  t,  f,g€  C G. 

tCG  \uCG  / 

Generally  since  G  is  not  necessarily  abelian,  fg  is  not  necessarily  equal  to  gf,  /,  g  G  CG. 
Denote  by  L(G;  C)  the  C-vector  space  of  all  C-valued  functions  on  G.  Every  /  €  L(G]  C) 
defines  a  formal  sum  in  C G,  /  =  Etec  /(*)*  and  we  can  identify  the  C-vector  space  CG 
with  the  C-vector  space  L(G;  C).  The  multiplication  fg  in  CG  corresponds  to  the  possibly 
noncommutative  convolution  in  L(G;  C) 

f*9(t)=  52  =  52  0(«)/(tu_1)»  />  5  €  L(G;  C). 

Under  this  identification  the  delta  function  8U>  u  e  G,  corresponds  to  a  formal  sum  having 
a  single  nonzero  coefficient  which  we  denote  by  u.  In  this  way  we  can  view  G  as  a  subset  of 
CG  and  a  basis  of  the  C-vector  space  CG. 

For  u  €  G  and  /  6  CG, 

«/  =  £/(*)«* 

tec  teG 

We  call  uf  the  left  translation  of  /  by  u.  Right  translation  can  also  be  defined  and  generally 
differs  from  left  translation,  unless  G  is  abelian.  For  u  e  G,  the  operator  L(u)  of  CG  defined 
by 

L{u)f  =  uf,  fe  CG, 

is  a  linear  isomorphism  of  the  C-vector  space  CG. 

A  subspace  V  of  the  C-vector  space  CG  is  called  G-invariant  if  for  all  u  G  G, 

UV  =  {«/  :  /  G  V}  C  F. 

A  G-invariant  subspace  V  of  CG  is  called  irreducible  if  the  only  G-invariant  subspaces 
of  V  are  (0)  and  V.  One  of  the  main  goals  of  nonabelian  group  harmonic  analysis  is  to 
characterize  the  G-invariant,  irreducible  G-invariant  and  direct  sum  decompositions  of  CG 
into  irreducible  G-invariant  subspaces.  For  an  abelian  group,  its  character  theory  provided 
all  the  necessary  tools  for  answering  these  questions.  Generally  these  problems  require  a  vast 
mathematical  machinery  for  their  solution.  However  for  the  nonabelian  groups  considered 
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in  this  text,  explicit  solutions  for  many  of  these  questions  will  be  derived  in  terms  of  abelian 
group  character  theory. 

For  g  E  CG,  the  operator  L(g)  of  C G  defined  by 

L(g)f  =  gf,  f  g  CG, 

is  a  linear  homomorphism  of  the  C-vector  space  CG.  Since 

L{g)  =  Y 9(t)L(t),  g  e  CG, 

tGG 

G-invariant  subspace  V  of  CG  are  CG-invariant  subspaces.  If  V  is  a  G-invariant  subspace 
of  CG  and  /  €  V,  then  for  all  g  E  CG,  gV  C  V. 

Idempotent  theory  provides  a  convenient  language  in  which  to  express  many  of  the  results 
in  the  following  chapters.  However  explicit  results  can  usually  be  written  in  terms  of  abelian 
group  character  concepts. 

A  nonzero  element  e  E  CG  is  called  an  idempotent  if  e2  =  e.  Two  idempotents  e\  and  e2 
in  CG  are  called  orthogonal  if  eie2  =  e2ei  =  0.  A  set  of  pairwise  orthogonal  idempotents 


{ej  :  1  <j<J} 


is  said  to  be  complete  if 


j 


3= 1 


Example  3  If  e  is  an  idempotent  in  CG,  then  {e,  1  -  e}  is  a  complete  set  of  orthogonal 
idempotents. 

Example  4  If  e  is  an  idempotent,  then  the  G-invariant  subspace  generated  by  e,  CGe,  has 
e  as  a  right  unit 

CGe  =  {a  E  CG  :  ae  =  o;}. 


Example  5  If  e  is  an  idempotent  then  CG  —  CGe  ©  CG(1  —  e). 

The  result  described  in  example  5  holds  generally  for  any  complete  set  of  orthogonal 
idempotents.  We  state  the  result  without  proof. 


Theorem  5  If 


{ej  :  1  <3  <J) 


is  a  complete  set  of  orthogonal  idempotents ,  then 


J 

CG  =  ^2  ©CGej. 
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Corollary  6  If  an  idempotent  e  can  be  written  as  the  sum 
e  =  ei  +  e2,  then 

C Ge  =  CGei  ®  C Ge2. 


of  two  orthogonal  idempotents, 


Example  6  If  CG  =  Wx  ©  W2,  where  Wi  and  W2  are  G-invariant  subspaces  and 

1  =  ei  +  e2 ,  e\  G  Wi,  e2  G  W2, 
then  {ei,e2}  is  a  complete  set  of  orthogonal  idempotents. 

Generally  we  have  the  following  result. 

Theorem  6  If  CG  is  the  direct  sum  of  G-invariant  subspaces,  CG  —  E/=i  ®Wvand, 

j 

1  =  ^2  ei >  eJ  ^  Wj  i  1  <  j  <  J, 

j= 1 

then  {e:  :  1  <  j  <  J}  is  a  complete  set  of  orthogonal  idempotents. 

Example  7  If  C Ge  —  W\  ®  W2,  the  direct  sum  of  G-invariant  subspaces  W\  and  W2  with  e 
an  idempotent  and  e  =  ei  +  e2,  ei  G  ,  e2  G  W2,  then  ei  and  e2  are  orthogonal  idempotents. 


Theorem  7  If  C Ge  is  the  direct  sum  of  G-invariant  subspaces 

j 

CGe  =  J2®WP 

3= 1 


with  e  an  idempotent  and 

e  =  J2 eJ>  e3  e  1  -  3  - 

3= 1 

then  {ej  :  1  <  j  <  J}  is  a  set  of  pairwise  orthogonal  idempotents. 

A  homomorphism  P  of  the  C-vector  space  CG  is  called  a  projection  if  P2  =  P.  Every 
subspace  W  of  CG  determines  a  projection. 

We  will  now  show  that  every  G-invariant  subspace  W  of  CG  is  generated  by  an  idem- 
potent.  Consider  any  projection  P  of  CG  satisfying  W  =  imP  and  define  the  mapping 
P0  :  CG  — ^  CG  by 

Po(ot)  =  Y2  u~1P{ua),  a  €  CG- 

u£G 

Po  is  a  homomorphism  of  the  C-vector  space  CG.  Since  P(ua)  G  W,  u  G  G  and  a  G  CG, 
and  W  is  G-invariant 


P0(a)  G  W,  a  G  CG. 
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Theorem  8  For  all  a  G  CG 
and  if  a  eW,  then 


Po(a)  =  aP0(l), 
P0(a)  =  a. 


Proof  For  t  G  G, 

t-'Po(t)  =  i  E  t-'u~'P(ut)  =  i  £  «-‘P(u)  =  P„(l), 

■'v  uGg  ueG 

by  a  change  of  variables.  The  linearity  of  Pq  proves  the  first  part.  If  ot  G  W ,  then  ucz  G  W , 
u£  G  and  P(ua)  =  ua  implying 

P0(ua)  =  Y,  u~lua  =  a» 

•'V  ixGG 

completing  the  proof. 

Since  im  P0  C  W  and  P0  acts  by  the  identity  mapping  on  W,  by  theorem  8  we  have  the 
following. 

Corollary  7  Pq  is  a  projection  of  CG  satisfying  W  =  ini  /  q  . 

Set  e  =  P0(l). 

Corollary  8  e  is  a  generating  idempotent  for  W , 

W  =  CGe,  e2  =  e. 


Since  every  G-invariant  subspace  W  has  a  generating  idempotent,  by  corollary  8,  there 
exists  a  G- invariant  subspace  W  such  that  CG  =  W  ©  W .  More  generally  we  have  the 
following  result. 

Corollary  9  If  Wi  and  W%  are  G-invariant  subspaces  of  CG  such  that  W\  C  W2  then  there 
exists  a  G-invariant  subspace  W[  ofW-2.  such  that  Wi  =  W\  ©  W[ . 

An  idempotent  e  G  CG  is  called  primitive  if  e  can  not  be  written  as  the  sum  of  orthogonal 
idempotents  in  CG. 

Theorem  9  e  is  a  primitive  idempotent  in  CG  if  and  only  if  CGe  is  irreducible. 

Proof  If  e  is  not  primitive  and  e  =  e\  +  e2,  where  e\  and  e2  are  orthogonal  idempotents, 
then  CG  =  CGei  ©  C Ge2  and  CGe  is  not  irreducible.  Conversely  if  IT  is  a  G-invariant 
subspace  of  CGe,  then  by  corollary  9,  CG  =  W  ©  W\  for  some  G-invariant  subspace  W'  of 
CGe.  Example  7  implies  e  is  not  a  primitive  idempotent,  completing  the  proof 

The  problem  of  constructing  G-invariant,  irreducible  G-invariant  and  direct  sum  decom¬ 
positions  of  CG  into  irreducible  G-invariant  subspaces  can  be  replaced  by  that  of  constructing 
idempotents,  primitive  idempotents  and  complete  sets  of  primitive  orthogonal  idempotents. 
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3.3.2  Group  algebras  of  A  B 

A  subgroup  A  of  a  group  G  is  called  normal  if  for  all  t  £  G, 

tAr1  =  {tar1  :  a  £  A}  C  A. 

G  is  said  to  be  the  semidirect  product  of  a  normal  group  A  and  a  subgroup  B  if  every  t£G 
can  be  written  uniquely  as 

t  —  xy ,  x  £  A,  y  £  B. 

In  this  case  we  write  G  =  A  <£  B. 

For  a  £  C A  and  (3  £  CB,  we  can  view  a,  {3  £  C G.  If  a/3  =  0,  then 

5Z  a(x)(3(y)xy  =  0 

x€A  y€B 

which  implies 

a(x)(3(y)  =  0,  x  £  A,  y  £  B, 

and  a  =  0  or  /3  =  0. 

Suppose  G  —  A  B  where  A  and  B  are  abelian  groups  of  orders  L  and  M.  Denote  the 
character  groups  of  A  and  B  over  C  by  A*  and  B*. 

B  acts  on  A*.  For  r  £  A*  and  y  £  B  define  ry  £  CG  by 

ry  =  yTy~x. 


Since  A  is  a  normal  subgroup  of  G 

ry  =  Yi,  rr{x)yxy~1  =  ^  r(y~1xy)x 

x&A  xGA 

is  also  in  A*. 

For  r  £  A* 

Bt  =  {tv  :  y  £  B} 

is  a  subset  of  A*  called  the  B-orbit  over  r.  A*  is  partitioned  into  the  disjoint  union  of  distinct 
B-orbits. 

Since  . 

L  r€,4*  M  A GB* 


we  have 


‘■jE  S>a, 


t£A *  X €B+ 

where  N  =  LM  is  the  order  of  G.  However  generally  r  and  A  do  not  commute.  We  do  have 
the  following. 

Theorem  10  For  t  =  xy£G,x£A,y£B, 

t  =  E  E  'r*{x-1)\(y-1)Tv\. 


t€A*  AGjB* 
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Proof  Since  xry  =  Ty{x~l)ry  and  yA  =  A (y  X)A, 

trX  -  xyrX  =  xryyX  =  Ty(x~1)X(y~1)TyX 


completing  the  proof. 

Since  G  is  a  basis  of  the  C-vector  space  C G,  by  theorem  10,  the  collection  of  products 
in  CG 

{rA  :  r  G  A\  X  G  B*} 

is  also  a  basis  of  CG.  Generally  CGtA  is  not  one-dimensional,  so  that  the  spaces 


{CGtA  :  r  G  A* ,  X  G  B*} 


intersect. 

Generally  to  form  direct  sum  decompositions  of  CG  into  left  G-invariant  subspaces  we 
must  modify  the  above  approach.  For  r  G  A* ,  define 

B(t)  =  {y  G  B  :  ry  =  r}. 

B(t)  is  a  subgroup  of  B  called  the  centralizer  of  r  in  B. 

Theorem  11  If  r 2  G  Bt\,  t\,  t2  G  A*,  then 

B(t1)  =  B(t2). 

We  can  assign  a  centralizer  in  B  to  every  B-orbit  in  A*.  The  assumption  that  B  is 
abelian  is  essential. 

Suppose  r  G  A*  and  consider  B(t).  For  y  G  B,  the  set 

yB(r)  =  {yz:z  G  B(r)}  C  B 

is  called  the  left  coset  of  B{r)  in  B  determined  by  y.  The  collection  B/B{t)  of  left  cosets  of 
B(t)  in  B  forms  a  partition  of  B.  A  set 

{ys  :  1  <  s  <  S}  C  B 

is  called  a  complete  system  of  representatives  of  B/B{r)  if  B  is  the  disjoint  union 

5=1 

The  B-orbit  Bt  has  order  S  and  Bt  =  {r^  :  1  <  s  <  S}.  Generally  the  y„  1  <  s  <  S 
depend  on  r  and  when  we  need  to  express  the  dependence  we  write  yj,  1  <  s  <  ST.  If  MT 
denotes  the  order  of  B(r), then  M  =  MTST,  where  M  is  the  order  of  B. 

We  will  now  show  that  the  collection  of  products 

~tA,  t  S  V.  A  €  B(t)\ 

lmt 

forms  a  complete  set  of  primitive  orthogonal  idempotents  for  CG.  We  break  up  the  proof 
into  the  following  three  theorems. 
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Theorem  12  For  r1;  r2  G  A*  and  Ax  G  A2  G  B(t2)* 

(tiAi)2  =  LMTit\X\ 


(r1A1)(r2A2)  =  0,  unless  tx  =  r2  and  Ax  =  A2 


Proof  Since  rf  =  Ti,  y  G  B(rx), 

nAi=  M2/)ri2/  =  I]  Ax(y)yri  =  Am. 

yeB(ri)  y€B(ri) 


Consequently 


(rxAx)2  =  Axt2Ax  =  LrxAf  =  LMTitx\x 


in  which  case 


completing  the  proof. 


(txAx)(t2A2)  =  Ai(r1r2)A2  =  0,  unless  rx  =  r2, 
(nAi)(riA2)  =  TiAiA2  =  0,  unless  Ax  =  A2, 


By  theorem  12,  the  collection  of  products  is  a  set  of  orthogonal  idempotents  in  C G. 


Theorem  13 


i  -  E  E 


Proof  The  theorem  follows  from 


t£A- 


X  =  iEr 

L  t<ea- 


1  =  T7-  E  tGAI*. 

1V1t  A€B(t)* 


Theorem  13  implies  completeness. 

For  r  G  A*  and  A  G  B(r)*,  by  theorem  12 


e  =  — tA 
LMt 


an  idempotent.  Since  j-r  is  an  idempotent  and  r A  =  At 


1  1 

e  =  yre  =  —  er. 

jL/  i-/ 
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Theorem  14  If  t  €  and  A  €  B(r)* ,  then  e  =  ^rA  is  a  primitive  idempotent. 

Proof  Assume  e  is  not  primitive  and  e  =  ei+e>2  where  t\  and  e 2  are  orthogonal  idempotents. 


Since 


ei  =  eie 


-eier 


-eir 


and 


we  have 


d  =  eei  =  —ree  1  =  —re\ 
L  L 


ex 


e? 


1  2 

^2relr  =  T^relr- 


I5' 


Since  A*  is  a  basis  of  CA,  we  can  write 


Consequently 

However 

and 

implying 


ci  =  tE  (  £  «i(^. vb'W  «i(Av)sc. 

L  yes  Vr'e^*  / 

ei  =  -^reir  =  ^^  £  ei(T'-  yjrr'yr. 

^  ^  yeB  r'eA* 

rr'yT"  =  0,  unless  r'  =  r  and  y  €  £(r), 

r2yr  =  L2ry,  y  e  B(t), 

=  JT  Y,  Ci(r,y)y  =  yrmb  mx  6  CS(r). 


ei 


yeB(r) 

The  same  argument  shows  that 


1 


implying 


e2  =  -TJTI2,  m2  €  C B(t), 


e  =  t^-tA  =  yr(mi  +  m2). 


LMT  L 

Since  t  G  CA  and  A,  m\  +  m2  €  C B(t), 


— -A  =  mi  +  m2. 
Mr 


We  will  show  that  mi  and  m2  are  orthogonal  idempotents  in  CB(r),  contradicting  the  fact 
that  -A- A  is  a  primitive  idempotent  in  CB(t). 


Since  rmi  =  mir, 


1  2  1 

-rmf  =  —rmi 
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implying  ml  —  m\.  The  same  argument  shows  m\  =  m2  and  mim.2  —  m2mi  —  0,  completing 
the  proof. 

We  have  proved  that  the  collection  of  products 

-i-rA  :  r  £  A* ,  A  £  B{t)* 

lmt 

is  a  complete  set  of  primitive  orthogonal  idempotents  for  C G  and  we  can  apply  the  results 
of  the  previous  chapter. 

Theorem  15 

c g=J2  J2  ®CGtX 

t€j4*  A €B(t)* 


with  the  direct  sum,  factors 


CGtX,  re  a*,  A  e  B{t)*, 

irreducible  G -invariant  subspaces  of  C G. 

If  G  is  not  abelian,  an  irreducible  G-invariant  subspace  of  C G  is  not  necessarily  one¬ 
dimensional.  For  groups  of  the  form  G  =  A  ^  B,  with  A  and  B  abelian,  we  will  determine 
bases  for  the  subspaces  CGrA,  r  £  A*,  A  £  B(r)*,  and  derive  fast  algorithms  for  relating 
components  in  the  G-basis  to  components  in  the  new  basis. 

Suppose  r  G  A*  and  A  €  B(r)*.  Set  ys  =  yTs,  1  <  s  <  S  =  ST  in  the  following  discussion. 
For  x  £  A  and  y  £  B,  with  y  =  ysz,  1  <  s  <  S  and  z  £  B(t), 


xyrX  =  xyszrX  =  XTy’yszX 
=  Ty’(x~1)X(z~1)ysTX. 

Theorem  16  For  r  £  A*  and  X  £  B(t)*,  the  set 

{ysT A  :  1  <  s  <  S} 

is  a  basis  of  the  C-vector  space  CGrA. 

Proof  Since  the  set 

{xyrX  :  x  £  A,  y  £  B} 


spans  CGrA,  the  set 

{ysrA  :  1  <  s  <  5} 

spans  CGrA.  We  must  show  that  this  set  is  linearly  independent.  Suppose  that 


0  =  ^a(s)ysr A,  a(s)  £  C,  1  <  s  <  S. 

3=1 
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Multiplying  on  the  left  by  ry‘,  1  <  t  <  S, 

s 

0  =  ifO  =  0i(s)TytTV3ysX 

5  =  1 

=  La(t)ytrX 

implying  a(t)  =  0.  Since  t  is  arbitrary,  1  <  t  <  S,  oc(t)  —  0,  for  all  1  <  t  <  S,  completing 
the  proof. 


For  a  E  C G, 


and  by  theorem  16 


L  tCA-  1VIt  AeB(r)* 


1  1 


sr 


TXTar A  =  E  ar\(s)yjTX,  aTx(s )  6  C. 
^  s~\ 


Theorem  17  For  a  €  C G,  r  €  A*  and  X  E  B(r)*, 


<*ta(s)  =  jjf  J2  (E  cx(xysz)Ty’(x  ^  A (z  1), 

L  Mt  zCB(t)  \xZA  ) 


where  ys  =  yTs ,  1  <  s  <  S . 

Proof  Since 

ST 

arX  =  E  E  E  a(xysz)xyszrX , 

5=1  zGB(t)  xGA 

the  theorem  follows  from 

xyszrX  =  Ty’ (x~l)X(z~1)ysTX>  z  E  B(t ),  1  <  s  <  ST- 


3.3.3  Convolution  Theorem 

Decompositions  of  KG  into  direct  sums  of  G-invariant  subspaces 

KG-Zew, 

3= 1 

lead  to  block  diagonal  matrix  representations  of  the  left  translation  operators  L(a),  a  E  KG, 
and  to  fast  algorithms  for  computing  products  in  KG. 

For  t  E  G,  the  matrix  P(t)  of  L(t)  relative  to  the  basis  G,  of  KG,  ordered  in  some  way,  is 
a  permutation  matrix  reflecting  the  group  structure  of  G.  The  matrix  S(a )  of  L(a),  a  E  KG, 
relative  to  G  is  a  linear  combination  over  K  of  the  permutation  matrices  P{t),  t  E  G.  If  G 
is  a  cyclic  group,  then  the  corresponding  matrices  P(t ),  t  E  G,  are  cyclic  shift  matrices  and 
the  corresponding  matrices  P(oi),  a  E  KG,  are  circulant  matrices.  In  this  case  the  direct 
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sum  decomposition  of  KG  given  by  the  characters  diagonalizes  the  circulant  matrices  with 
the  Fourier  transform  describing  the  change  of  basis. 

Generally  since  Wj  is  G-invariant,  L(a),  a  £  KG,  maps  Wj  into  itself.  Denoting  by 
Tj(a)  the  matrix  of  the  restriction  of  1(a)  to  some  basis  of  Wj,  the  matrix  direct  sum 

rW  =  EJ)(“). 

3= 1 

is  the  matrix  representation  of  L(a)  relative  to  the  basis  of  KG  formed  by  piecing  together 
the  bases  from  the  Wj,  1  <  j  <  J.  We  say  the  new  basis  is  compatible  with  the  direct  sum 
decomposition  of  KG. 

To  compute 

otp  =  L(a)0,  a,  0  £  KG, 

we  write  (3  in  terms  of  the  compatible  basis,  compute  T(a)  and  form  the  matrix  product  of 
T(a)  with  the  coordinates  of  (3  in  the  compatible  basis  and  then  translate  the  result  back  to 
the  basis  G. 

In  the  following  sections,  we  will  carry  out  the  program  for  the  direct  sum  decomposition 

KG=Y,  E  ®KGt\, 

A€B(t)* 


relative  to  the  idempotent  basis 

yTsrX,  t  £  A*,  X  £  B(t)*,  1  <  s  <  St. 

Matrix  Representation  of  L 

For  t  £  A*  and  A  £  B(r)*,  denote  by  LTX(g),  g  £  KG,  the  restriction  of  L(g)  to  KGtX  and 
by  TrX(g)  the  matrix  of  Lr\(g)  relative  to  the  idempotent  basis  of  KGtX.  The  matrix  of 
L(h)  relative  to  the  idempotent  basis  neglecting  order)  is  given  by  the  matrix  direct  sum 

E  E  Tr»(«). 

7 €A*  AeB(r)* 


Right  translation  from  A 

Suppose  t  £  A*  and  A  £  B(t)*.  Set  ys  =  yTs,  1  <  s  <  S  =  ST.  For  x  £  A, 


LT\{x)y)sTX  = 


xysrX 

XTV3ysX 

Tys(x~1)Ty“ysT 

Ty’{x~l)ysTX, 


1  <s<S, 


proving  the  following  result. 
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Theorem  18  For  r  E  A*, 

TrX(x )  =  diag(Ty,(x~1))i<s<s,hpp\  E  x  E  A, 

with  ys  =  yTs,  1  <  s  <  S  =  ST. 

The  diagonal  matrix  in  the  theorem  depends  solely  on  r  and  is  independent  of  A  leading 
to  the  following  corollary. 

Corollary  10  For  r  E  A*, 

Y  ®diag  (ry‘  (x"1 ))  ,  x  E  A, 

A€B(T)*®rr*(z)=/MT 

with  ys  =  yTs,  1  <  s  <  S  =  ST. 

The  corollary  describes  the  matrix  of  the  restriction  of  L(x),  x  E  A,  to  the  G-invariant 
subspaces  J2xe b(t)*  ®FCGt A. 

Corollary  11 

T(x)  =  Y  ®  dia9  (rW'(x_1))i  s  ,  x  E  A, 

TGA* 

where  rVs  =  ty*}  1  <  s  <  S  =  ST. 


Left  translation  from  B 

The  description  of  T(t),  t  E  B,  requires  deeper  analysis.  For  r  E  A*  and  A  E  B{r)*,  we 
describe  the  matrix  TrX(t),  teB,hy  the  following  steps. 

•  TrX(yTs),  1  <  s  <  5r. 

.  TrX  (yrsB(T)) ,  1  <  s  <  ST. 


Since  B  is  the  disjoint  union 

B  =  Ui  <s<SryTsB{r) 

all  TrX(t),  te  B,  are  described.  The  steps  are  dependent  on  r  E  A*. 
Set  ys  —  yj,  1  <  s  <  ST.  For  1  <  t  <  S, 


ytys  =  y'sz>  l  <  s,  s'  <  S,  z  e  B(t). 

There  exists  a  permutation  nj  E  Permm(S)  such  that 

ytVs  =  y^{s)zt(s),  zt(s )  E  B(t),  1  <  s  <  S. 

For  1  <  t  <  S, 

ytysr\  =  y^(s)Zt{s)rX 
=  y^(s)TZt{s)X 

=  X(zt(s)~1)y7rj (s)xX,  1  <  sS. 

Denote  by  P[  the  5x5  permutation  matrix  defined  by 

p;  =  [e„.W  ■  ■  ■  «..(«)]  .  1  <  <  <  S- 
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Theorem  19  For  t  £  A*  and  A  G  B(t)*, 

TTX(.yn  =  P;diag( 

The  permutation  7rtr  and  the  permutation  matrix  P(T,  1 <  t  <  ST,  depend  solely  on  r  6  A* 
and  are  independent  of  A  €  B(r)*,  leading  to  the  following  corollary. 


Corollary  12  For  r  €  A*, 

£  ©TrA(2/tT)  =  (/MT®PtT)  £  ^(A((^(s))"1)1<s<sr' 

A€B(t)*  A€B(t)* 

By  the  corollary,  the  matrix  of  the  restriction  of  L(yJ)  to  the  G- invariant  subspace 

£  0/sTGrA 

A  €B(r)* 

uses  the  same  permutation  matrix  PtT  on  each  KGtX,  A  G  B(r)*.  This  step  requires  for  each 
t  €  A*,  ST  permutation  matrices  P[,  1  <  t  <  ST. 

Suppose  there  exists  y  =  yT  G  B  such  that 

ys,  0  <  s  <  S  =  ST 


is  a  complete  system  of  left  coset  representatives  of  B/B (r).  The  idempotent  basis  of  KGtX 


is  given  by 


ysr A,  0  <  s  <  S. 


For  0  <  t  <  S, 

ytysT  A  =  yt+srA, 

where  t  +  s  is  taken  modulo  S.  We  now  have 


TM)  =  P\ 


where  P  is  the  S-point  cyclic  shift  matrix 

f  1  0  1 


L  o  10 

For  z  G  B(t), 


zysrX  =  A [z  ')y stX , 

1  <  s  <  5, 

implying  that  Tt\(z)  is  the  scalar  matrix 

Trx{z )  -  A(z-1)7St, 

z  G  B(t). 
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Theorem  20  For  r  G  A*  and  A  €  B(r)*, 

Trx(ytz)  =  A(z-1)PtTcfeap  , 

with  1  <  t  <  ST  and  z  G  B(t). 

Corollary  13 

E  ©7V»(»z)  =  (■?«,»  ^7)  E  A(z-1)rfia9(A((z,(s))-1)1Si£Sr. 

AGB(t)*  A€B(t)* 

3.4  Characters  of  Nonabelian  Groups 

Assume  throughout  that  the  characteristic  of  K  is  relatively  prime  to  the  order  |G|  of  G. 
3.4.1  Basic  Properties 

A  mapping  p  :  G  — >  Kx  is  called  a  one- dimensional  character  of  G  over  K  if  p  is  a  group 
homomorphism, 

p(xy)  =  p(x)p{y),  whenever  x,y  eG.  (1) 

One-dimensional  characters  of  nonabelian  groups  share  with  abelian  group  characters  the 
property  that  they  generate  one-dimensional  G-invariant  subspaces.  Since  we  will  not  use 
the  more  general  concept  of  a  character  of  a  group  representation  of  G,  we  will  use  the  term 
character  to  mean  one-dimensional  character. 

As  a  function  on  G,  a  character  p  of  G  can  be  viewed  as  an  element  in  KG, 

p(x)x ■ 

x£G 

Several  of  the  proofs  below  are  exactly  the  same  as  in  the  abelian  case,  but  we  include  them 
for  completeness. 

Theorem  21  Suppose  p  is  a  character  of  G.  For  any  t  E  G  and  a  €  KG} 


tP= pt  =  p(t  x)p 


and 


where 


ap  =  pa  =  a(p)p, 

&{p)  = 
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Proof  Prom  ( 1)  and  changing  variables,  we  have 

tp  =  Yj  p(x)tx  =  X  p(t~lx)x  =  p(rl)  X  P(x)x  =  p^~x)p' 

xeG  xgg  xgg 

and 

pt  =  X  Pfa)^  =  X  pO^-1)®  =  p(t_1)  X  P(x)x  =  p(*  )p> 

x€G  x€G 

proving  the  first  statement. 

Using  the  first  statement,  we  have 

ocp  =  X  a(*)*p  =  (x  a(*)p(r1))  />> 

teG  VteG  / 

and  ,  \ 

pa  =  Y,  <*$)&  =  (  X  a(*)P(rl) )  P» 

teG  VteG  / 

completing  the  proof. 

For  7  g  KG,  define  the  centralizer  G( 7)  of  7  in  G  by 

G(7)  =  {teG:tj  =  7 1}. 

(5(7)  is  a  subgroup  of  G.  The  center  of  KG  is  the  set  of  all  7  £  K G  such  that  G  G(7)- 
7  is  in  the  center  of  KG  if  and  only  if  we  have  70  =  07,  for  all  a  e  XG.  By  Theorem  10, 
every  character  p  of  G  is  in  the  center  of  KG. 

Characters  of  G  generate  one-dimensional  G-invariant  subspaces.  For  a  character  p  of  G, 

KGp  =  {S(p)p  :  a  €  KG}.  (2) 


In  fact  we  have  the  following. 

Theorem  22  Every  one- dimensional  G-invariant  subspace  of  KG  is  generated  by  a  char- 
acter  of  G  over  K . 

Proof  Suppose  KGj,  7  €  KG,  is  one-dimensional.  For  t  £  G, 

tj  =  aft) 7,  aft)  e  Kx. 

Without  loss  of  generality,  we  can  assume  7(1)  =  1.  Since  (st) 7  =  sft'y),  we  have  a(st)  = 
a(s)aft),  for  all  s,teG.  a  and  7  are  then  characters  of  G  over  K,  completing  the  proof. 

Denote  by  G*  the  collection  of  all  characters  of  G  over  K . 
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Theorem  23  For  characters  pi  and  pi  of  G  over  K,  if  pi  ^  pi,  then 


p\  =  \G\Pi 


and 


P1P2  —  P2P1  —  0. 


Proof  Since  p\  =  pi(pi)pi,  we  have  by  ( 1), 

PM  =  Hpi(*)pi(t_1)  =  LM1)  =  lGl> 

tGG  tGG 

the  first  result  follows. 

If  Pi  ¥"  P2,  P2{pi)  =  Pi(P2)  and 

P1P2  =  P2i.P1)  Pi  =  Pi{p2)p2 

imply  pip2  =  0,  completing  the  proof. 

As  a  consequence  of  the  theorem,  if  pi  and  p2  are  distinct  characters  of  G  over  K, 

Pi  (Pa)  =  Y,Pi{t)p2(t~1)  =  0. 

tGG 


By  Theorem  13,  we  have  the  following. 

Corollary  14  The  collection 

is  a  system  of  primitive  pairwise  orthogonal  idempotents  in  KG. 

The  corollary  implies  G*  is  linearly  independent. 

For  a  finite  nonabelian  group  G,  KG  cannot  decompose  into  the  direct  sum  of  one¬ 
dimensional  G-invariant  subspaces,  since  this  would  imply  that  G  is  isomorphic  to  a  group 
of  diagonal  matrices  which  is  clearly  an  abelian  group.  Unlike  the  finite  abelian  group  case 
splitting  over  K ,  the  collection  in  Corollary  (14)  is  not  complete.  However  a  complete  system 
of  orthogonal  primitive  idempotents  in  KG  exists.  Necessarily  some  of  these  idempotents  are 
not  characters  of  G  and  generate  G-invariant  subspaces  of  dimension  greater  than  one.  In 
the  following  section,  for  the  classes  of  groups  considered,  we  will  compute  such  systems.  We 
will  lay  the  foundation  for  doing  so  by  first  establishing  additional  properties  of  characters 
of  G  and  characters  of  subgroups  of  G. 

For  y  e  G  and  a  €  KG,  define  ay  e  KG  by 

-1  =  a{y~lxy)x. 


ay  =  yay 


(3) 
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The  mapping  of  KG  defined  by  ( 3)  is  an  automorphism  of  the  group  algebra  KG.  By 
Theorem  10,  if  p  is  a  character  of  G  then  for  all  y  E  G,  py  —  p. 

A  subgroup  H  of  G  is  called  normal  if  for  all  y  E  G, 

yHy~l  =  {yxy~l  :  x  E  H}  C  H. 

If  a  G  KH, 

a  =  Y  a(x)x 

xCH 

and  y  E  G  then  since  yxy~l  E  H,  for  all  a :  E  H, 

ocy  =  Y  0'(x)yxy~1  =  Y  Oi{y~1xy)x 

xCH  x£H 

is  again  in  KH.  The  automorphism  (3)  maps  KH  onto  itself,  for  every  y  E  G. 

Suppose  throughout  that  H  is  a  normal  subgroup  of  G.  For  a  character  p  of  H  we  have, 
for  every  y  G  G,  that  py  is  a  character  of  H,  i.e.,  the  mapping 

x — »  p(y~1xy),  x  G  H, 

is  a  character  of  H.  In  this  way  G  acts  on  H*,  the  collection  of  characters  of  H  over  K. 
This  action  plays  a  major  role  throughout  this  work.  We  will  introduce  certain  constructions 
associated  to  this  action  at  this  time  and  leave  to  later  sections  their  application. 

Consider  p  G  H*  and  define 

Gp  =  {pv  :  y  G  G}. 

Gp  is  a  subset  of  H*  called  the  G-orbit  over  p.  Two  G-orbits  which  intersect  must  be  equal. 
As  a  result,  H*  is  the  disjoint  union  of  its  distinct  G-orbits.  A  subset  of  H*, 

{pi,  •  •  •  ,Pr} 

is  called  a  complete  set  of  representatives  for  the  G-orbits  in  H*  if  H*  is  the  disjoint  union 

H*  =  uf=1G/,r. 

For  p  G  H*,  the  centralizer  of  p  G  G, 

G(p)  —  {y  E  G  :  py  =  p} 


is  a  subgroup  of  G  containing  H . 

Theorem  24  For  p  G  H*  and  y  G  G,  we  have 

G(py)^yG(p)y-1. 
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Proof  If  z  6  G(p)  then  since 

(: yzy~l)pv  =  yzy~lypy~l  =  yzpy~l  =  ypzy~l 
=  ypy-'yzy-1  =  /?( yzy -1), 

we  have  yzy~l  G  G(py).  A  similar  argument  shows  that  u  G  G{py)  implies  y~xuy  €  G{p), 
completing  the  proof. 

In  general,  G(p)  is  not  a  normal  subgroup  of  G.  However,  if  it  is,  we  have  the  following 
corollary. 

Corollary  15  If  p  G  H*  and  G(p)  is  a  normal  subgroup  of  G,  then  for  all  z  G  G,  G(pz)  = 
G(p). 

Some  more  group  theory  is  necessary  for  the  next  set  of  results.  For  y  G  G,  the  subset  of 

G, 

Hy  =  {zy  :  z  g  H) 

is  called  a  left  H -coset  over  y.  Left  H- cosets  which  intersect  are  equal  and  so  G  is  the  disjoint 
union  of  its  distinct  left  H- cosets.  Since  H  is  a  normal  subgroup  of  G,  a  group  multiplication 
can  be  placed  on  the  collection  of  left  H- cosets  by 

(Hyi)(Hy2)  =  H(yiy2 ),  Vi,  2/2  €  G. 

Normality  is  used  to  show  that  the  multiplication  is  well  defined  in  the  sense  that  if  Hyi  — 
Hy[  and  Hy2  =  Hy'2  then  H{yxy2)  =  H(y[y'2),  y[,  y'2  G  G.  The  resulting  group  is  denoted 
by  H\G.  H  is  the  identity  element  and  Hy~x  is  the  inverse  of  Hy  in  H\G. 

A  subset  of  G 

{yT  ■  0  <  r  <  -R},  (4) 

is  called  a  complete  set  of  representatives  for  H\G  if  G  is  the  disjoint  union  of  the  left-cosets 
over  yr,  0  <  r  <  R.  We  will  always  take  yo  =  1  in  G. 

3.4.2  Primitive  Idempotents 

A  character  p  of  II  over  K  determines  the  primitive  idempotent  jj^p  in  KH.  However,  in 
general,  jj^p  is  not  a  primitive  idempotent  in  KG.  We  continue  assuming  that  H  is  a  normal 
subgroup  of  G. 

Theorem  25  Suppose  p  G  H*  and  set  f  =  If 

f  —  e  l  +  e2i 

where  e\  and  e2  are  orthogonal  idempotents  in  KG,  then  e\  and  e2  are  in  the  KG(p). 
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Proof  Since  ex  and  e2  are  orthogonal  idempotents, 

fe  i  =  e\  +  e2ei  =  ex  =  e\  +  eie2  =  ex/. 


From  ( 4).  we  can  write 

R- 1 

ej  =  Y,  Y,  ej  (xyr)xyr. 

r=0  zeH 

Theorem*  10  implies 

j  fi-i 

ex  =  /ex  =  tttt  E  Z  ei(a:?/r)(pa;)r/r 
I"  I  r=0 

=  /  53  f  J2  eifayrXz-1))  Vr  =  /  5Z  ai(02/r,  e  K. 

r=0  \x£H  /  ’’=0 

Since  „  . 

1  fi-i 

ex  =  ex/  =  777^  a(0PPVryr 

1^1  r=0 

we  have  from  Theorem  13,  that 

ex  =  f  <x(r)yr, 

yrQG(p) 

which  since  H  C  G(p)  implies  ex  €  KG(p).  The  same  argument  shows  that  e2  G  KG(p), 
completing  the  proof. 

Normality  of  H  in  G  is  essential  since  it  is  required  for  ip-py"1  €  H*  and 

PVrPVr1  =  0 

unless  yT  G  G(p). 

Corollary  16  If  p  €  H*  and  H  =  G(p)  then  j j^p  is  a  primitive  idempotent  in  KG. 

Theorem  25  and  Corollary  16  motivate  a  strategy  for  constructing  primitive  idempotents 
in  KG.  For  the  groups  considered  in  the  following  two  sections  this  strategy  will  lead  to 
algorithms  for  constructing  complete  system  of  primitive  pairwise  orthogonal  idempotents  in 
KG.  We  begin  with  a  character  p  of  H  over  K  and  search  for  extensions  of  p  to  characters 
Pi  over  subgroups  Hx  in  G  containing  H  with  the  hope  that  Hx  =  G(px). 

Assume  throughout  that  H  is  a  normal  subgroup  of  G  such  that  the  group  of  left  H -cosets 
H\G  is  an  abelian  group.  Under  this  assumption  any  subgroup  Hx  of  G  containing  H  is  a 
normal  subgroup  of  G. 

Suppose  Hi  is  a  subgroup  of  G  containing  H.  Since  H  is  the  identity  element  in  the 
group  H\G ,  for  all  xx,  x2  €  G,  the  commutator 

x^x^x^1  G  H. 
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If  x2  E  Hi,  then 

xix2x{x  E  Hx 2  C  Hi,  X\  E  G, 
proving  Hi  is  a  normal  subgroup  of  G. 

Suppose  p  is  a  character  of  H  over  K.  A  subgroup  Hi  of  G  admits  an  extension  of  p  if 
H  C  Hi  and  there  exists  a  character  pi  of  Hi  over  K  extending  p, 

pi(x)  =  p(x),  x  E  H. 

We  say  simply  that  pi  is  an  extension  of  p  to  Hi. 

Theorem  26  If  Hi  admits  an  extension  of  p  then  Hi  C  G(p) . 

Proof  Suppose  pi  is  an  extension  of  p  to  Hi.  For  x  E  H  and  y  E  Hi, 

Pi(yx)  =  Pi{y)p{x) 


and  since  yxy~x  E  H, 

pi{yx)  =  pi{yxy~ly)  =  p{yxy~x)pi{y)  =  py{x)pi{y). 

As  a  result 

p(x)  =  py(x),  x  E  H,  y  E  Hu 

proving  the  theorem. 


Theorem  27  Suppose  Hi  admits  an  extension  of  p.  If  pi  and  p2  are  two  extensions  of  p  to 
Hi,  then 

Hi  c  G{pi)  =  G(P2)  C  G(p). 


Proof  For  y  E  G(pi)  and  c  E  Hi, 

Pi(yxy~lx~l)  =  pX(x)pi{x~x)  =  1. 


Since  yxy  lx  1  E  H, 

p(yxy~1x~1)  =  1. 

Writing 

P2(yxy~x)  =  p2{yxy~lx~lx)=  p{yxy~xx~l)p2{x) 

=  P2{x), 

we  have  G(pi)  C  G(p2).  The  same  argument  shows  G(p2)  C  G'(pi),  proving  G(pi)  —  G(p2). 
If  we  take  x  E  H,  then  yxy~x  and  x~x  are  in  H  and 

p(yxy~xx~x)  =  py(x)p(x~x),  x  E  H, 
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implying  G(pi)  C  G(p)  and  completing  the  proof. 

Suppose  Hi  admits  an  extension  p\  of  p.  Denote  by  p{Hx)  the  collection  of  all  extensions 
of  p  to  H\.  We  can  construct  p(Hi)  from  pi  and  the  collection  of  all  characters  of 

the  finite  abelian  group  H\Hi  over  K. 

For  x  €  (H\Hi)*,  the  mapping  x  '■  Hi  — »  K ,  defined  by 

x(®)  =  x(Hx),  X  e  Hi, 

is  an  extension  of  the  trivial  character  on  H  to  Hi, 

x{x)  =  1,  for  all  x  6  H. 

We  can  identify  ( H\Hi )*  with  the  collection  of  all  extensions  of  the  trivial  character  on  H 
to  Hi. 

Suppose  pi  is  an  extension  of  p  to  Hx.  For  any  x  €  {H\HX)*,  the  mapping  pf  :  Hx  — > 
Kx,  defined  by 

Pi(x)  =  x(x)Pi(x)>  x  e  Hi, 

is  an  extension  of  p  to  Hx.  Conversely  if  p2  is  any  extension  of  p  to  Hx  then 

x(x)  =  Pi{x)P2ix~1),  x  e  Hi, 

is  an  extension  of  the  trivial  character  on  H  to  H\,  proving  the  following  result. 

Theorem  28  Suppose  pi  is  an  extension  of  p  to  Hi.  Then 

p(Hi)  =  {p?  :  X  €  (H\Hi)*}. 

Every  extension  of  p  to  Hx  has  a  unique  representation  of  the  form 

pf,  xe(H\Hi)*. 


In  particular,  the  order  of  p{Hx)  is  the  order  of  {H\HX)*.  If  H\HX  splits  over  K,  the 
order  is  the  order  of  H\HX. 

We  will  use  Theorem  28  to  prove  a  result  which  is  essential  for  establishing  completeness 
of  certain  systems  of  primitive  pairwise  orthogonal  idempotents. 

Theorem  29  Suppose  Hx  admits  an  extension  of  p  and  H\HX  splits  over  K .  Then 

1  _  1 

W\p-mei£j„f" 
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Proof  Suppose  p[  is  an  extension  of  p  to  Hv  By  Theorem  28,  for  any  x  6  Hu 

£  piW  =  p;w  £  *(#*)■ 

Piep(HO  xew\H ,)- 

The  right  hand  summation  vanishes  unless  x  G  H  and  we  have  for  x  E  H, 

Pi€p(H,)  1  1 

completing  the  proof. 

Suppose  Hi  admits  an  extension  px  of  p.  Choose  a  complete  set  of  representatives 

{yT  :  0  <  r  <  R} 

for  the  group  H\H\.  Every  y  €  Hi  has  a  unique  representation  of  the  form  y  =  xyr ,  x  6  H, 

0  <  r  <  R.  We  can  write 

K— 1  H-l 

Pi  =  Y  Y  Pl{xyr)xyr  =  Y  P(X)X  Y  Pl(Vr)yr 

r=0  x£H  x€H  r=0 

and  we  have 

R — 1 

Pl  =  P  Y  Pl(3h)yr- 

r= 0 

In  particular,  if  p 1  is  a  character  of  H  over  K  and  p*  ^  p,  then  (J p\  =  p\pf  =  0.  This  result 
has  the  following  important  consequence. 

Theorem  30  For  distinct  characters  p  and  p'  of  H  over  K  and  extensions  pi  and  p[  to  Hi 
and  H[  we  have  pip\  =  P1P1  =  0. 

By  Theorem  26  and  the  corollary  to  Theorem  ??,  if  H  =  G{p)  then  p  cannot  be  extended 
and  |^|P  is  a  primitive  idempotent  in  KG.  Explicit  algorithm  for  constructing  primitive 
idempotents  is  presented  in  Subsection  4.1.1. 

3.5  A  Theorem  of  Mackey  and  Wigner 

Suppose  G  is  a  finite  group  having  the  form  G  =  A  <$.  H  where  A  is  an  abelian  group  of 
order  |Aj  splitting  over  the  field  K  and  H  is  an  arbitrary  subgroup.  The  main  result  in 
this  section  is  a  theorem  due  to  Wigner  and  Mackey.  The  result  provides  a  framework  for 
constructing  complete  systems  of  primitive,  pairwise  orthogonal  idempotents  for  semidirect 
product  groups.  The  construction  of  such  systems  in  the  preceding  section  can  be  viewed  as 
a  special,  simple  case. 

The  Mackey-Wigner  (M-W)  little  group  theorem  is  not  the  end  of  the  story  even  for 
semidirect  product  groups  since  it  requires  knowledge  of  complete  systems  of  primitive, 
pairwise  orthogonal  idempotents  for  certain  subgroups  of  H.  In  the  preceding  section,  H  is 
abelian  and  such  systems  are  given  in  terms  of  abelian  group  character  theory.  We  will  use 
the  M-W  theorem,  along  with  the  results  of  section  3.4,  to  find  such  systems  for  KG,  where 
H  has  the  form  H  =  B  $.C,  with  B  and  C  abelian  groups. 
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3.5.1  M-W  Theorem 

Suppose  for  r  G  A*, 

{/rT  :  1  <  r  <  rT} 

is  a  complete  system  of  primitive,  pairwise  orthogonal  idempotents  for  KH(r). 

Then  the  set  of  products 

iV/;’  TeA’’l-r5r~  (5) 

is  a  complete  system  of  primitive,  pairwise  orthogonal  idempotents  for  KG. 

Theorem  31  Suppose  r  €  A*  and  a  G  KG  satisfies 

Then  a  =  ^nn,  m  G  KH{r). 

Proof  Write 

«=EE  a(xy)xy. 

y£H  x£A 

Since  a  =  ^ ra,  we  have 

a  =  TT?  E  (E  a(zy)Tz)  V  =  |4r  E  (E  (x(xy)T(y-1))  y, 

\A\  yCh  J  VCH  \xCA  / 

which  we  write  as 

Q  =  tt\t  E  ®(v)y>  5(y) G  K 

K1!  yCH 

Since  a  =  a  (  j^t) ,  we  have 

a  =  ni2  E  «(y)rryy- 

I"!  yCH 

ttv  =  0,  unless  y  G  #(r),  implies 

a  =  TTTr  E  «(y)y=nTTm>  rn&KH(r), 

\A\  yCH(r)  I"' 

completing  the  proof. 

In  the  following,  we  will  several  times  use  the  result  that  if  a  G  KA  and  (3  e  KH,  the 
a(3  =  0,  if  and  only  if  a  —  0  or  (3  =  0. 

Theorem  32  Suppose  r  G  A*  and  f  is  a  primitive  idempotent  in  KH(t).  Then  e  =  rW 
is  a  primitive  idempotent  in  KG. 
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Proof  r  and  /  commute,  implying  e  is  an  idempotent.  Suppose  e  is  not  primitive  and 
we  can  write  e  =  e\  +  e2,  where  e\  and  e2  are  orthogonal  idempotents  in  KG.  Since 
jLre  -  e  =  e.  (pfj'r)  and  eei  =  t\  =  eie,  we  have 


mTei = ei = 61 


Applying  Theorem  31, 


61  =  Rirmi’ 


mi  G  KH(t). 


In  the  same  way  e2  =  ,  m2  G  KH{r).  Since  r  and  mi  commute, 


ei 


|A|Tm 


2 

i) 


1  1 

By  the  observation  preceding  the  theorem,  mi  =  m\.  Arguing  in  the  same  way,  mi,  m2  are 
orthogonal  idempotents  in  KH(r).  Since 

iV/=iV(mi+m2)’ 

the  same  observation  implies  /  =  mi  +  m2,  a  contradiction,  proving  the  theorem. 


3.5.2  Proof  of  the  M-W  theorem 

By  Theorem  32,  \X\Tfr  >  t  G  A*,  1  <  r  <  rr,  is  a  primitive  idempotent  in  KG.  Since  r  and 
fj  commute  the  collection  of  products  ( 5)  is  a  system  of  primitive  orthogonal  idempotents 
for  KG.  Since 


MI&S*  14  .S-T 


(£/; 

\r=l 


=  m.S-t  =  1' 


the  system  is  complete,  proving  the  theorem. 


3.6  Abelian  Extension 

The  M-W  theorem  will  be  used  to  construct  an  algorithm  for  computing  complete  system 
of  primitive  orthogonal  idempotents  for  abelian  extensions.  We  begin  by  assuming  that  A 
is  a  subgroup  of  a  semidirect  product  G  =  A  <£  B  of  finite  abelian  groups  A  and  B.  By  the 
M-W  algorithm,  we  can  build  a  complete  system  of  primitive  orthogonal  idempotents  for  G 
from  the  characters  of  A.  The  approach  can  be  tailored  to  the  subgroup  A. 


AF98T009  Phase  II  Final  Report:  Mathematical  Results 


40 


The  projection  P  :  G  — >  B  defined  by 

Px  =  z ,  xeG,x  =  yz,yeA,zeB, 

is  a  group  homomorphism  of  G  onto  B  having  kernel  A.  Since  A  C  4  ^  PA,  where  PA  is 
a  subgroup  of  B,  we  can  assume  without  loss  of  generality  that  B  =  PA. 

Suppose  p  is  a  character  of  A.  Define 

{Pi  =  {z  G  B  :  [z,  A]  C  ker(p)} 

and 

Ai  =  {x  G  A  :  Px  G  Pi}. 

Pi  is  a  subgroup  of  P  and  Aj  is  a  subgroup  of  A  containing  A  H  A. 

For  X\  G  Ai,  x  =  yz,  y  G  A  and  z  G  Pi,  define 

pi(*)  =  p(v)- 

Theorem  33  px  is  a  character  of  Ay  extending  the  restriction  of  p  to  A  fl  A. 

Proof  Suppose  xx,  x2  G  Ai  and  write  x\  =  yxZx,  x2  =  y2^2,  Vi,  V2  €  A,  zu  z2  G  Pi-  Since 

xxx2  =  [zi,xi]yiy2z\Z2 

and  \z\ ,  x2]  G  ker(p), 

Pi{x  ix2)  =  pi(xi)pi(x2), 

proving  pi  is  a  character  of  Ai.  Clearly,  pi  extends  the  restriction  of  p  to  A  D  A. 

The  critical  condition  used  in  this  theorem  is  that 

Ai  cA$Bu 

with  [P1}  Ai]  C  ker(p)  which  implies  that  pi  is  a  character  of  Ai.  In  general  pi  is  not  a 
maximal  extension  in  A  of  the  restriction  of  p  to  A  n  A  For  this  to  occur  we  require  that 

Ai  =  {x  G  A  :  [x,  Ax]  C  ker(p)}. 

Define 

P2  =  {z  G  P  :  [z,  Ax]  C  ker(p)} 

and 

A2  =  {x  G  A  :  Px  G  P2}. 

P<2  is  a  subgroup  of  P  containing  Pi  and  A2  is  a  subgroup  of  A  containing  Ax- 
Theorem  34 


A2  =  {x  G  A  :  [x,  Ax]  C  ker(p)}. 
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Proof  Suppose  x  G  A2  with  x  =  yz,  y  e  A,  z  e  P2.  For  any  *1  G  Ax,  with  *1  =  yi2i> 
yi  G  A,  zi  G  Bu  we  have 

[x,  zi]  G  ker(p),  [z,  xx]  G  fcer(p) 


implying 

[x,*i]  =  [z,xi][x,zi]  G  fcer(p). 

Since  xj  is  arbitrary  in  Ax,  [x,  Ai]  G  ker(p)  and  we  have  shown  that 

A2  C  {x  G  A  :  [x,  Ai]  C  ker(p)}. 

Conversely,  take  x  G  A  such  that  [x,  Ai]  C  ker(p).  For  any  xi  G  Ai,  [x,xx]  G  ker(p). 
Since  [x,zi]  G  ker(p),  and 

[x,xi]  =  [z,xi][x,zi]  G  ker(p), 

we  have  [z,Xi]  G  ker(p)  and  z  G  B2.  This  shows  that  x  G  A2j  completing  the  proof. 


If  Ai  =  A2,  then  px  is  a  maximal  extension  in  A  of  the  restriction  of  p  to  A  D  A  and 
jA-jpj  is  an  idempotent  in  CA.  Otherwise,  define 

B3  =  {z  G  B  :  [z,  A2]  C  ker(p)} 

and 

A3  =  {x  G  A  :  Px  G  B3}. 

Since 

[£i,A2]  C  [Pi,  A]  C  ker(p), 
we  have  that  Bx  is  a  subgroup  of  B3.  In  the  same  way, 

[£s,  Ax]  C  [P3,A2]  C  ker(p), 
implies  B3  is  a  subgroup  of  P2.  We  have  shown 

Ai  c  A3  c  A2  C  A. 

Replacing  A  by  A2  and  Ai  by  A3  we  can  continue  in  this  way  and  construct  a  subgroup  BM 
of  B  such  that  the  subgroup  Am  of  A  defined  by 

Am  =  {*£  £  A  :  Px  G  Pm} 

satisfies  the  two  conditions 

[Am,  Pm]  C  ker{p) 

Am  =  {x  G  A  :  [x,  Am]  C  ker(p)}. 

The  first  implies  that  the  mapping  pm  ■  A m  — >  Cx  defined  by 

Pm(x)  =  p{y),  x  G  Am,  x  =  yz,  y  G  A,  z  G  BM 

is  a  character  of  Am  extending  the  restriction  of  p  to  A  0  A.  The  second  shows  that  pM 
is  a  maximal  character  extension  in  A  of  the  restriction  of  p  to  A  D  A  and  is  an 

idempotent  in  CA. 
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3.7  Equivalence 

Many  of  the  most  important  results  on  the  structure  of  the  algebra  CG  are  described  in 
terms  of  two-sided  G- invariant  subspaces  of  CG.  We  will  see  that  CG  decomposes  into  the 
direct  sum  of  all  its  irreducible  two-sided  G-invariant  subspaces.  This  is  the  first  step  in 
proving  the  Wedderburn  structure  theorem  which  constructs  a  faithful  representation  of  CG 
as  the  direct  sum  of  complete  matrix  algebras  over  C. 

The  decomposition  of  CG  into  the  direct  sum  of  its  irreducible  two-sided  G-invariant 
subspaces  is  necessarily  unique.  The  analogous  statement  for  irreducible  G-invariant  sub- 
spaces  does  not  hold,  but  will  up  to  equivalence  between  such  spaces.  Equivalent  irreducible 
G-invariant  subspaces  are  the  building  blocks  for  the  irreducible  two-sided  G-invariant  sub¬ 
spaces. 

These  results  extend  the  classical  DSP  algorithm  for  computing  cyclic  convolutions  by 
fast  Fourier  transforms  and  diagonal  matrix  multiplications.  For  nonabehan  groups,  the 
diagonal  matrix  multiplications  are  replaced  by  block-diagonal  matrix  multiplications.  This 
result  is  called  the  Wedderburn  Structure  theorem  which  has  many  generalizations. 

For  a  e  CG,  right-translation  by  a  is  the  C-linear  transformation  R(a)  of  CG  defined 

by 

R(oc)f  =  fa,  f  €  CG. 

By  the  associative  law  in  CG,  R(a)  commutes  withe  all  left- translations  of  CG, 

R{a)Q3f)  =  /3R{a)f,  0,  f  e  CG.  (6) 

In  fact  any  mapping  T  of  CG  into  itself  satisfying  (6)  is  right-translation  by  T(l), 

T(f)  =  T(f  •  1)  =  fT(  1),  /  e  CG. 

Suppose  W  is  a  G-invariant  subspace  of  CG.  For  a  6  CG,  denote  by  Rw{&)  the 
restriction  of  R(af  to  W .  R\y(a)  is  a  C-linear  transformation  of  IF  onto  IFo;, 

Wol  =  {fa  :  /  €  W}. 

Wot  is  a  G-invariant  subspace  of  CG  by  (6).  It  may  be  the  {0}  subspace. 

Suppose  e  is  a  generating  idempotent  for  the  G-invariant  subspace  W.  A  mapping  T  of 
W  into  CG  satisfying 

T(0f)  =  0T{f),  P,  /  e  W, 

equals  Rw(oi)  where  a  =  T(e). 

The  following  results  of  Schur  describe  right-translations  on  irreducible  G-invariant  sub¬ 
spaces  and  serve  as  important  tools  throughout  this  chapter. 

Theorem  35  If  W  is  an  irreducible  G-invariant  subspace  and  a  €  CG,  then  Wa  —  {0} 
or  Wa  is  an  irreducible  G-invariant  subspace.  IfWa  ±  {0},  then  Rw{ot )  is  a  C-linear 
isomorphism  ofW  onto  Wa. 
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Proof  Suppose  Wa  ±  {0}  and  Rw(a)  is  not  a  C-linear  isomorphism.  Then 

V  =  {/  g  W  :  fa  =  0}  ±  {0} 

is  a  G-invariant  subspace  of  W.  By  irreducibility,  V  =  W  and  Wo  —  {0},  a  contradiction, 
proving  Rw( a)  is  a  C-linear  isomorphism  of  W  onto  Wa. 

Suppose  Wa  ±  {0}  and  U  is  a  G-invariant  subspace  of  Wa.  Then 

{feW-.faeU} 

is  a  G-invariant  subspace  of  W .  It  is  either  the  {0}  subspace  and  U  =  {0}  or  it  is  W  and 
JJ  =  Wa,  proving  W a  is  irreducible  and  completing  the  proof. 


Corollary  17  IfWa  ±  {0}  and  e'  is  a  generating  idempotent  for  Wa,  then  there  exists  a 
unique  h  EW  such  that  e'  =  ha.  We  have 

•  Rw«(h)Rw(a)  =  Iw, 

•  Rw{a)Rw<x{h)  =  Iwcn 

•  (Wa)h  =  W. 


Proof  Since  Rw(&)  is  a  C-linear  isomorphism  of  W  onto  Wa,  there  exists  a  unique  heW 
satisfying  e'  =  ha.  To  complete  the  proof  we  need  only  to  prove  the  second  statement.  Take 
g  €  Wa.  Then  ge'  —  g  and 

(gh)a  =  g{ha)  =  ge'  =  g, 

completing  the  proof. 

For  an  irreducible  G-invariant  subspace  W  of  CG  we  can  have  W  =  Wa  for  some 
a  €  CG.  For  example,  if  e  is  a  generating  idempotent  for  W,  then  W  =  We  and  Rw(e)  —  Iw- 
We  have  the  following  result  of  Schur. 

Theorem  36  Suppose  W  is  an  irreducible  G-invariant  subspace  of  CG.  If  a  €  CG  such 
that  W  =  Wa,  then  there  exists  a  €  Cx  such  that 

Rw(a)  =  alw. 


Proof  By  Theorem  35,  Rw{a)  is  a  C-linear  automorphism  of  W.  Since  C  is  algebraically 
closed,  Rw{a)  has  an  eigenvalue  a  £  Cx  and  there  exists  a  nonzero  /0  €  W  such  that 

foa  =  a/o- 


Consider 

W{a)  =  {/eW:  Rw(a)f  =  a/}. 
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For  (3  G  CG  and  f  eW(a),  since 

(Pf)a  =  (3{fa )  =  atff), 

we  have  0f  e  W(a)  and  W(a)  is  a  G-invariant  subspace  of  W.  Since  W{a)  +  {0},  W{a)  = 
W ,  completing  the  proof. 

Corollary  18  If  e  is  a  generating  idempotent  for  W ,  then  for  some  a  €  C  , 

eoc  =  eae  =  ae. 


Suppose  W  is  an  irreducible  G-invariant  subspace  of  CG.  A  G-invariant  subspace  W'  f 
{0}  is  called  equivalent  to  W  if  there  exists  a  e  CG  such  that  W'  =  W a.  By  Theorem  35 
and  its  corollary,  if  W'  is  equivalent  to  W,  then  W'  is  irreducible  and  W  is  equivalent  to 
W' .  In  this  way  we  define  an  equivalence  relation  on  the  collection  of  irreducible  G-invariant 
subspaces  of  CG.  Denote  by  WR  the  collection  of  all  irreducible  G-invariant  subspaces 
equivalent  to  W . 

We  have  the  following  criterion  for  equivalence. 

Theorem  37  Suppose  W  and  W'  are  irreducible  G-invariant  subspaces  having  generating 
idempotents  e  and  e' .  Then  W  and  W'  are  equivalent  if  and  only  if  there  exits  a  nonzero 
7  G  CG  satisfying 

7  =  eq  =  je1  =  eje1 . 

In  this  case  W'  =  Wj. 

Proof  Suppose  W  and  W'  are  equivalent  and  W'  —  Wot,  a  G  CG.  Since  W  =  We  and 
Wot  =  ( Wa)e we  have  Wa  =  Weae'.  Setting  7  =  eae', 

7  =  ej  =  7e'  =  eje1  7^  0, 

and  W'  =  Wj. 

Conversely  if  there  exists  7  7^  0  in  CG  satisfying 

7  =  e7  =  7e'  =  e'ye ', 

then  W7  C  W .  W 7  is  a  nonzero  subspace  of  W'  and  W'  =  Wy,  completing  the  proof. 

We  need  assume  in  the  theorem  that  7  =  eye1  ±  0  since  7  =  eye'  implies  7  =  e7  and 
'"y  =  7e;. 

For  the  following  discussion,  we  write  U  ~  V  to  mean  that  the  irreducible  G-invariant 
subspaces  U  and  V  are  equivalent.  Consider  a  direct  sum  decomposition  of  a  G-mvanant 
subspace  W  of  CG, 

W  =  £©Wt, 

t=i 

into  irreducible  G-invariant  subspaces  Wu  1  <  t  <  T.  We  will  show  that  up  to  equivalence 
such  decompositions  are  unique.  The  main  result  needed  is  the  following. 
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Theorem  38  If  V  is  an  irreducible  G -invariant  subspace  contained  in  W ,  then 

VC  ©Wt, 

wt~v 

the  direct  sum  taken  over  all  Wt  equivalent  to  V ,  1  <  t  <  T. 

Proof  Suppose  e  is  a  generating  idempotent  for  V  and  write 

r 

e  =  Y^at,  OLtC  Wt. 

t= i 


We  have  that 

vcYl  ©w. 

OLt  ^0 


Since  e  is  an  idempotent 


Uniqueness  of  representation 


T 

e  =  e2  =  ^  eat,  eat  €  Wt. 

t=i 

in  the  direct  sum  implies 


eat  =  at,  1  <  t  <  T. 


For  any  1  <  t  <  T  such  that  at  ^  0,  we  have  eat  /  0  and 

{0}  ±  Uat  C  Wt. 


Irreducibility  implies 


Vat  =  Wt, 


and  Wt  ~  V,  completing  the  proof. 


at  7^0 


Corollary  19  Suppose  Vs,  1  <  s  <  5,  is  a  set  of  equivalent  irreducible  G -invariant  subspaces 
contained  in  W .  Then  5Zf=i  Vs  C  Yiwt~Vi  ©Wt. 

Consider  any  set  of  irreducible  G-invariant  subspaces  VSl,  1  <  Si  <  *?i  such  that  every 
Vs,l<s<S,  is  equivalent  to  a  unique  V^,  1  <  Si  <  S\. 


EKc  E  ©wt. 

v.~vai  Wt~vtl 


Corollary  20  For  1  <  Si  <  Si, 
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By  the  proceeding  theorem  and  its  corollaries,  if  W  has  two  direct  sum  decompositions 
into  irreducible  G-invariant  subspaces 

T  S 

W  =  j>Wt  =  £©14, 

t- 1  s= 1 

up  to  rearrangement  of  the  factors,  W.,  is  equivalent  to  Vs,  1  <  s  <  S  and  by  dimension 
S  <T.  Reversing  the  roles,  we  have  the  following  result. 

Theorem  39  Suppose  W  =  £?=i  Wt  =  £f=i  ©14  are  direct  sum  decompositions  of  a  G- 
invariant  subspace  W  into  irreducible  G-invariant  subspaces  Wt,  1  <  t  <  T,  and  Vs,  1  <  s  < 
S.  Then  T  —  S  and  up  to  rearrangement,  Wt  is  equivalent  to  Vt,  1  <  t  <  T. 

3.7.1  Two-Sided 

A  G-invariant  subspace  J  of  CG  is  called  two-sided  if  J  is  invariant  under  all  right-translations, 

Jot  C  <7,  ot  E  CG. 

A  nonzero  two-sided  G-invariant  subspace  J  is  called  irreducible  if  the  only  two-sided  G- 
invariant  subspaces  of  J  are  {0}  and  J.  An  arbitrary  two-sided  G-invariant  subspace  par¬ 
titions  the  collection  of  all  irreducible  two-sided  G-invariant  subspaces  into  two  sets;  those 
contained  in  J  and  those  having  trivial  intersection  with  J . 

Suppose  W  is  a  G-invariant  subspace  and  J  and  J'  are  two-sided  G-invariant  subspaces. 
WJ  is  a  two-sided  G-invariant  subspace  contained  in  J  and  JW  is  a  G-invariant  subspace 
contained  in  J  n  W.  In  particular  J  J'  is  two-sided  and  is  contained  in  J  l~!  J  .  For  a  direct 
sum  J  ©  J'  we  have  J  J'  —  J'J  C  J  fl  J'  =  {0} . 

Suppose  W  is  an  irreducible  G-invariant  subspace  of  CG  and  J  is  a  two-sided  G-invariant 
subspace.  WR  denotes  the  collection  of  all  irreducible  G-invariant  subspaces  equivalent  to 
W.  For  a  €  J,  we  have  either  Wa  =  {0}  or  Wa  G  WR.  Since  Wot  G  J, 

W J  C  U w'€Wr(J  Fl  VF  )  (7) 

In  particular  if  J  fl  W'  =  {0}for  all  W'  E  WR,  then  W  J  =  {0}. 

Theorem  40  Suppose  J  is  a  two-sided  G-invariant  subspace  and  W  is  an  irreducible  G- 
invariant  subspace  such  that  J  GW  ^  {0} .  Then 

W'  C  J,  for  all  W'  E  WR. 


Proof  Since  J  GW  is  a  G-invariant  subspace  of  W  and  J  GW  ^  {0},  we  have  JGW  -W 
and  W  C  J.  The  invariance  of  J  under  right-translations  implies  Wa  C  J,  for  all  a  E  CG, 
completing  the  proof. 
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From  the  theorem  a  two-sided  G-invariant  subspace  J  partitions  the  set  of  all  equivalence 
classes  of  irreducible  G-invariant  subspaces  into  two  subsets,  those  having  all  its  elements 
contained  in  J  and  those  having  all  its  elements  having  trivial  intersection  {0}  with  J. 

For  an  irreducible  G-invariant  subspace  W  of  CG,  define 

J(w)=  £ 

WZWR 

J(W)  is  G-invariant  and  is  the  minimal  G-invariant  subspace  containing  all  W'  G  WR. 

Theorem  41  J(W)  is  an  irreducible  two-sided  G-invariant  subspace. 

For  a  G  CG  and  W'  G  WR ,  we  have  W'a  =  {0}  or  W'a  G  WR.  From 

J(W)a=  Y,  W'aC  £  W'  =  J(W),  ol  G  CG, 

W'ewR  w'ewR 

J(W)  is  a  two-sided  G-invariant  subspace. 

Suppose  J  ^  {0}  is  a  two-sided  G-invariant  subspace  contained  in  JiW).  If  e  is  a 
generating  idempotent  for  J{W), 

J  =  JecJJ(W)  ±  {0}. 

Assume  J  Pi  W'  =  {0},  for  all  W'  G  WR.  Then 

JW'  C  J  n  W'  =  {0},  for  all  W'  G 
and  JJ(W)  =  {0},  a  contradiction.  We  must  have 

J  n  W'  7^  {0},  for  some  W'  G  WR, 

which  by  Theorem  40  implies  J  =  J(W).  J{W)  is  irreducible,  completing  the  proof. 

By  Theorem  39,  J{W )  has  an  essentially  unique  direct  sum  decomposition  into  irreducible 
G-invariant  subspaces 

J(W)  =  Y®V°-  (8) 

5—1 

Choose  T  <  S  such  that,  perhaps  up  to  reordering  of  the  factors  in  (8),  every  Vs,  1  <  s  <  S, 
is  equivalent  to  exactly  one  Vu  1  <  t  <  T.  We  will  show  T  =  1.  Denote  by  Jt  the  direct 
sum  of  all  the  factors  of  (8)  equivalent  to  Vu  1  <  t  <  T.  Then 

T 

J(W)  =  £®Jt- 

1 

We  will  show  that  J\  is  a  two-sided  G-invariant  subspace  which  by  irreducibility  of  J{W) 
implies  J\  —  J(W)  and  T  —  1. 
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Theorem  42  J{W)  is  the  direct  sum  of  irreducible  G -invariant  subspaces  all  of  which  are 
equivalent  to  W . 

Proof  Suppose  V  is  any  factor  in  the  given  direct  sum  decomposition  of  Jx .  Since  J (W)a  C 
J(W),  for  all  a  €  CG, 

T 

Va  C  for  all  a  G  CG. 

t=i 

By  Theorem  38,  Va  =  {0}  or  Va  is  equivalent  to  Vi,  implying 

Va  C  Ju  for  alia  G  CG. 

Since  V  is  an  arbitrary  factor  of  the  direct  sum  decomposition  of  Ji, 

J\Oc  C  Ji,  for  all  a  G  CG, 

and  Ji  is  a  two-sided  G-invariant  subspace,  completing  the  proof. 


Corollary  21  Every  irreducible  G-invariant  subspace  of  J{W)  is  equivalent  to  W. 
Corollary  22  Suppose 

{Ws  :  1  <  s  <  S} 

is  a  set  of  pairwise  inequivalent  irreducible  G-invariant  subspaces.  Then  the  sum  £f=1  J{WS) 
is  a  direct  sum 

Proof  By  the  preceding  theorem,  for 

j(ws)J(Wt)  c  J(WS)  n  J(Wt)  =  {0}. 

If 

0,  OLt  C  J(Wt), 

t=l 

and  es  is  an  idempotent  generator  for  J(WS),  1  <  s  <  S,  then 


ases  =  as  =  0, 


completing  the  proof. 


Suppose  J  is  a  two-sided  G-invariant  subspace  and  E(  J)  denotes  the  collection  of  all 
irreducible  two-sided  G-invariant  subspaces  contained  in  J.  By  Theorem  40,  there  exists  a 
set 

{Wt :  1  <  t  <  T} 

of  irreducible  G-invariant  subspaces  contained  in  J  such  that  every  irreducible  G-invariant 
subspace  contained  in  J  is  equivalent  to  exactly  one  Wt,  1  <  t  <  T,  and  J(Wt )  C  J, 
1  <  t  <  T. 


AF98T009  Phase  II  Final  Report:  Mathematical  Results 


49 


Theorem  43  J  =  Y%= 1  ©«/(Wi)  where 

E(J)  =  {J{Wt)  :  1  <  t  <  T}. 


Proof  By  Corollary  18  of  Theorem  42, 

t=  i 

Suppose  J'  is  an  irreducible  two-sided  G-invariant  subspace  contained  in  J.  If  J  n  J(Wt)  = 
{0},  for  all  1  <t<T,  then  JJ'  =  {0},  a  contradiction.' We  must  have 

J'  n  J(Wt)  ±  {0},  for  some  1  <  t  <  T. 

Irreducibility  implies  J'  =  J(Wt),  completing  the  proof. 


Corollary  23  J  is  the  direct  sum  of  all  the  irreducible  two-sided  G-invariant  subspaces 
contained  in  J. 

Applying  Theorem  43  and  its  corollary  to  CG  we  have  the  following  result  which  we  state 
as  a  theorem  because  of  its  importance. 

Theorem  44  CG  is  the  direct  sum  of  all  its  irreducible  two-sided  G-invariant  subspaces. 

Corollary  24  If  J  is  a  two-sided  G-invariant  subspace  of  CG,  then  J' ,  the  direct  sum  of 
all  the  irreducible  two-sided  G-invariant  subspaces  having  trivial  intersection  {0}  with  J,  is 
the  unique  two-sided  G-invariant  subspace  of  CG  satisfying  CG  =  J  ©  J  ■ 

Corollary  25  A  two-sided  G-invariant  subspace  J  ±  {0}  has  a  unique  idempotent  generator 
e  and  e  is  central  in  CG. 

Proof  By  the  preceding  corollary,  CG  =  J  ©  J' .  Write 

1  =  e  +  e',  eeJ.e'e  J' . 

e  and  e'  are  orthogonal  idempotents  generating  J  and  J'  as  G-invariant  subspaces.  For 
g  e  CG, 

g  =  g  ■  1  =  ge  +  ge!  =  1  •  g  =  eg  +  e'g. 

Since  J  and  J'  are  two-sided,  ge,  eg  G  J  and  ge' ,  e'g  G  J' ■  Uniqueness  of  representation 
implies  eg  =  ge  and  e  is  central  in  CG. 

Suppose  e\  is  also  an  idempotent  generator  of  J.  Since  e  is  central, 


e  =  eei  =  eie  —  e\ 
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completing  the  proof. 

Suppose  J  is  a  two-sided  G-invariant  subspace  having  a  direct  sum  decomposition 

J  =  X>VS, 

S=1 

into  irreducible  G-invariant  subspaces  Vs,  1  <  s  <  S.  As  in  the  discussion  leading  up 
to  Theorem  42,  after  rearrangement  if  necessary,  there  exists  T  <  S,  such  that  every  Vs, 
1  <  s  <  S,  is  equivalent  to  exactly  one  V,  1  <t<T,  and 

J  =  f>J(Vt). 

t= l 

Denoting  by  Jt  the  direct  sum  of  all  the  factors  Va,  1  <  s  <  S,  equivalent  to  Vu  we  have 

T 

t= l 

In  particular  Jt  is  an  irreducible  two-sided  G-invariant  subspace. 
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4  Development  of  Fast  Algorithms 

4.1  Character  extension 

4.1.1  Construction  of  Primitive  Idempotents 

Assume  that  H  is  a  normal  subgroup  of  G  and  H\G  is  an  abelian  group  splitting  over  the  field 
K.  Suppose  p  is  a  character  of  H  over  K.  By  Theorem  16  and  the  corollary  to  Theorem  25, 
if  H  =  G(p)  then  p  cannot  be  extended  and  j j^p  is  a  primitive  idempotent  in  KG.  Suppose 
H  is  strictly  contained  in  G(p).  The  following  algorithm  constructs  a  nontrivial  extension  of 

P- 

H  is  strictly  contained  in  G(p). 

•  Choose  y  G  G(p)  and  y^H  and  form  H i  the  subgroup  generated  by  H  and  y. 

•  Find  the  smallest  positive  integer  R  such  that  yR  G  H. 

•  Compute  A  =  p{yR)- 

•  Find  any  R- th  root  a  of  A  in  K,  aR  =  A,  possible  since  H\G  splits  over  K. 

•  Write  every  x  G  Hi  uniquely  in  the  form 

x  =  uyr ,  uGff,  0  <  r  <  R. 

•  Construct  the  mapping  pi  :  Hi  — >  K'x  by 

pi(uyr)  =  p{u)ar ,  0  <  r  <  R. 


Pi  is  an  extension  of  p  to  H\. 

By  repeated  use  of  the  algorithm,  possible  since  every  subgroup  of  G  containing  H  is 
normal  in  G,  we  can  construct  a  subgroup  M  of  G{p)  admitting  an  extension  7  of  p  to 
M  such  that  M  -  G( 7).  In  this  case  ^7  can  be  extended  no  further  and  is  a  primitive 
idempotent  in  KG. 

Theorem  28  can  then  be  used  construct  the  collection  p(M)  of  all  extensions  of  p  to  M. 
•  For  each  x  €  ( H\M )*,  form  the  mapping  7*  :  M  — >  Kx,  by 

7x(x)  =  x(ifx)7(x),  x  G  M. 

The  collection  p(M )  of  all  extensions  of  p  to  M  is  given  by 

p(jM)  =  {7X  :  X  ^  (H\M)*}. 

By  Theorem  23  and  the  corollary  to  Theorem  ??,  the  collection 

{jsp :  7  6 
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is  a  system  of  primitive,  pairwise  orthogonal  idempotents  in  KG  satisfying,  by  Theorem  29, 

—  1  V 

For  each  character  p  of  H  over  K,  denote  by  Mp  any  subgroup  of  G{p)  which  admits  an 
extension  q  of  p  satisfying  Mp  —  G ( 7 ) .  Denote  by  p(Mp)  the  collection  of  all  extensions  of 
p  to  Mp.  Then 

I(p)  =  :  7  e 

is  a  system  of  primitive,  pairwise  orthogonal  idempotents  for  KG  satisfying 

W\P==  ^ 

\H\  fel(p) 


By  Theorem  30, 

I  =  Upe  H'I(P) 

is  a  system  of  primitive,  pairwise  orthogonal  idempotents  for  KG  satisfying 

W\  ^  p  = 

1 73  I  /€/ 


For  H  an  abelian  group, 

1  =  t4[  E  P 

l"l  p€H* 

and  the  above  construction  leads  to  a  complete  system  of  primitive,  pairwise  orthogonal 
idempotents  for  KG. 

Up  to  scalar  multiple,  the  primitive  idempotents  of  KG  resulting  from  the  algorithm 
are  characters  of  subgroups  of  G.  If  H  is  abelian,  then  a  complete  system  of  primitive, 
pairwise  orthogonal  idempotents  of  KG  can,  up  to  constant  multiple,  be  given  by  characters 
of  subgroups  of  G.  If  H  is  not  abelian  this  need  no  longer  be  the  case. 


4.1.2  Coefficient  computation 
For  y  €  B,  define  ay  €  C  A  by 

Qy(x)  =  a(xy),  x  e  A, 

and  compute  the  Fourier  transform  of  ay  over  A, 

«3/(r)  =  7  Y,  Oiy(x)r{x-1),  t  e  A*. 

M  Fourier  transforms  over  A  are  required  at  this  stage.  Place  the  results  into  a  two- 
dimensional  array  over  A*  x  B 

6t(r,y)  =  ay(r),  r  G  A* ,  y  €  B. 
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For  each  r  G  A*  we  will  compute  the  coefficients 

«ta(s),  1  <  s  <  ST,  A  G  B(t)*. 

For  1  <  s  <  ST,  define  7J  €  CB(r)  by 

7l(*)  =  7 {Ty,,ysz),  z  G  -B(r),  ys  =  yj, 

and  compute  the  Fourier  transform  of  7J  over  B(r) 

■i:w  =  TT  E  Aesw. 

z6B(t) 

5t  Fourier  transforms  over  B(t)  are  required  in  this  stage.  Implementing  this  stage  as  r 
runs  over  A*  requires  EreA*  ST  Fourier  transforms  over  abelian  groups  over  varying  sizes. 
By  theorem  17,  for  each  r  G  A* , 

arA(s)=7j(A),  A  G  B(t)*,  1  <  s  <  St. 


First  Stage  a  G  C  G 

oty  G  C A,  y  E.  B 

FT  over  A,  y  G  B. 

&t  G  CA*,  y  G  B 

a  G  L(A*  x  J5) 

Second  Stager  G  >1* 

a  G  L(A*  x  B) 

7 J  G  CB(r),  1  <  s  <  5t 

FT  over  J5(r),  1  <  s  <  5T. 

7J  G  CB(t)*,  1  <  s  <  5t 


<»«(«)  =  71(A),  AB(r)*.  1  <  »  <  Sr. 
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4.2  Algorithm  G  =  A  <$  (B  <$  C) 

For  r  €  A* 

•  Compute  by  the  primitive  idempotent  construction  algorithm  a  complete  system  of 
primitive,  pairwise  orthogonal  idempotents  for  KH(t), 

fr,  1  <r  <rT. 

•  Form  the  products 

r/;,  1  <  r  <  rT. 


The  collection  of  products 

jijT r  €  A-,1  <  r  <  rT, 

is  a  complete  system  of  primitive,  pairwise  orthogonal  idempotents  for  KG. 

Up  to  a  constant  multiple,  the  idempotents  of  KH{t),  t  €  A*,  constructed  in  the  al¬ 
gorithm  are  characters  of  H(t).  As  a  result,  up  to  a  constant  multiple,  the  system  of 
idempotents  of  KG  resulting  from  the  algorithm  consists  of  characters  of  subgroups  of  G 
containing  A.  This  is  the  case,  even  though  A\G  is  not  abelian. 

By  the  Mackey- Wigner  (M-W)  little  group  theorem,  we  can  accomplish  the  first  task  by 
designing  an  algorithm  which,  for  each  r  G  A*,  computes  a  complete  primitive  idempotent 
system  in  KH(r).  In  the  next  section,  we  develop  an  algorithm  for  this  purpose  based  on 
the  algorithms  for  A  §B.  As  we  will  see,  combined  with  the  M-W  algorithm,  this  leads 
to  primitive  idempotents  of  KG  which  are  characters  of  subgroups  of  A  •$.  H(t).  We  make 
critical  use  of  this  property  in  the  design  of  the  basis  and  expansion  coefficient  algorithms 
established  in  section  4. 

4.2.1  Idempotent  Systems  for  Subgroups  of  H  =  B  <£  C 

Suppose  A  is  a  subgroup  of  H .  Since  B  and  C  are  abelian  groups,  A  n  B  is  an  abelian, 
normal  subgroup  of  A  and  A / A  Pi  B  is  an  abelian  group.  The  algorithms  for  previous  efforts 
apply. 

Algorithm  I  Complete  primitive  idempotent  system  in  A" A 
A 

For  p  e  (A  Pi  B)* 

•  Construct  a  subgroup  Ap  of  A  admitting  an  extention  7  of  p  such  that  ^7  is  a 
primitive  idempotent  in  K A. 

•  Construct  extpAp,  the  collection  of  all  extensions  of  p  to  Ap. 

Idem(p)  =  ^  e  extpAp}. 


•  Form  the  set 
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Idem(p)  is  an  idempotent  system  in  AA  satisfying 


1 

|Bn  A|p 


1 


E  e. 

eGldem(p) 


B 


•  Implement  A  for  all  p  G  (B  D  A)*. 

•  Form  the  set 

UpG(BnA)*  Idem(p). 

Since  A  Pi  B  is  an  abelian  group  splitting  over  A ,  the  collection 

Upe(snA)*  I  dem(p) 

is  a  complete  primitive  idempotent  system  in  AT  A. 

Combining  this  algorithm  with  the  M-W  algorithm,  we  have  the  following  algorithm  for 
constructing  an  idempotent  system  for  AG. 


Algorithm  II  Complete  primitive  idempotent  system  in  AG 
A 

For  r  e  A*, 

•  Apply  algorithm  I  to  construct  a  complete  primitive  idempotent  system  in  KH{r). 

fr,  1  <r  <rT. 

•  Compute  the  products 

rfr,  1  <r  <rT. 

B 

•  Implement  A  for  all  r  €  A*. 

•  Form  the  set 

{rfr  :  r  €  A*,  1  <  r  <  rT}. 

By  the  M-W  algorithm  this  set  is  a  complete  primitive  idempotent  system  in  AG. 

To  establish  notation  for  the  following  sections,  we  will  explicitly  describe  A  in  algorithm 
II. 

Algorithm  III  Step  A  of  algorithm  II 
A’ 

For  pe  (H(t)  fl  B)*, 

•  Construct  a  subgroup  H{t)p  of  Hir)  admitting  an  extension  7  of  p  such  that 
is  a  primitive  idempotent  in  KH{t). 
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•  Construct  the  collection  ext pH{t)p  of  all  extensions  of  p  to  H{r)p. 

•  Form  the  set  .. 

IdemT(p)  =  { |  ,  r— 1<$  :  S  £  ext pH{r)p}. 

\hKt)p\ 

B’ 

•  Implement  A’  for  all  p  €  (H(r)  H  B)* . 

•  Form 

Idem(r )  =  yJpe(H{r)nByIdemr{p). 

Idem(r)  is  a  complete  primitive  idempotent  system  in  KH(t). 

4.2.2  Basis  Algorithm 

The  primitive  idempotents  in  AG  resulting  from  algorithm  II  are  constant  multiples  of 
characters  of  subgroups  of  G.  The  basis  algorithm  developed  in  this  section  and  the  expansion 
coefficient  algorithm  of  the  next  section  significantly  depend  on  this  form  of  the  primitive 
idempotents. 

Suppose  t  £  A*,  p£  (B  fl  H(t))*  and  H(r)p  and  ext pH(t)p  are  constructed  by  algorithm 
III.  We  will  design  an  algorithm  for  constructing  a  basis  for  the  irreducible  G-invariant 
subspaces  KGr'y,  7  G  ext pH{r)p. 

Choose  a  complete  system  of  representatives 

ui,  1  <  Z  <  L, 

for  the  space  of  right  H{r)p- cosets  in  H.  Every  y  £  H  has  a  unique  representation  of  the 
form 

y  =  Uiw,  1  <  l  <  L,  w  £  H{r)p. 

In  general,  v-i,  1  ^  l  ^  L,  and  L  depend  on  t  and  p  but  for  the  most  part  we  suppress  this 
dependence  in  the  following  discussion. 

For  t  £  G,  we  can  uniquely  write 

t  =  xy  —  xuiw,  x  £  A,  y  £  H,  1  <  l  <  L,  w  £  H{t)p. 

tT  7  =  XUl'UJT'Y  =  XUiTUf^UiTW  7 

=  Ti(x~l)j(w~i)uiTJ, 

where  77  =  1  <  l  <  L.  Since  the  set  {tr-y  :  t  £  G}  spans  KGtj,  the  set  {u(T7  :  1  < 

l  <  L]  spans  KGt 7. 

Theorem  45  Suppose  t  £  A* ,  p  £  ( B  D  H(r ))  and  U[,  1  <  Z  <  T,  is  o  complete  system  of 
representatives  for  the  H{r)p-cosets  in  H.  Then 

{u;T7  :  1  <  l  <  L} 

is  a  basis  of  KGrj,  for  every  7  G  ext pH(r)p. 
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Proof  We  must  prove  linear  independence.  Suppose 

L 

^a(Z)uiT7  =  0,  a(l)  €  K. 

i= i 

For  1  <  m  <  L, 

L 

Tm  X  oL{i)nua  =  o, 

i=i 

where  n  =  uitu[1  ,  1  <  l  <  L.  Since  rmr,  =  0  unless  rm  =  r/,  in  which  case  =  |A|rm,  we 
have 

X)  a{l)uiTj  =  0, 

where  the  sum  is  over  all  1  <  l  <  L  such  that  €  #(t).  We  have  used  the  equivalence 

between  rm  =  r,  and  u^uiTuf1^  =  r.  Multiplying  by  u~\ 

u"1  X  a{l)uiTj  =  X  a(0ru'1u/7  =  0. 

V‘Tn1'U'l€:H(T)  Urr}  UiGH(t) 

Since  H(t)p  is  a  normal  subgroup  of  H(t)  and  is  the  centralizer  of  7  in  H{t),  we  have, 
arguing  as  above,  that 

X  ol{1)uit 7  =  0. 

This  sum  has  only  one  term  implying 

a(m)umTo/  =  0 

and  a(m)  =  0.  Since  1  <  m  <  L  is  arbitrary,  a(m )  =  0,  1  <  m  <  L,  completing  the  proof. 

In  the  following  algorithm,  we  suppress  the  construction  of  centralizer  subgroups. 

Algorithm  IV  Basis  for  KGt A 
For  r  £  A*  and  p  £  (B  Pi 

•  Construct  H(t)p  and  ext pH(t)p  by  algorithm  III. 

•  Construct  U(,  1  <  l  <  L,  a  complete  system  of  right  H(t)p- cosets  in  H. 

For  7  £  ext pH(t)p> 

•  Compute  the  products 

UlTJ,  1  <  l  <  L. 

By  the  theorem,  {itjry  :  1  <  l  <  L}  is  a  basis  of  KGtj. 
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4.2.3  Expansion  Coefficient  Algorithm 

Suppose  r  £  A*,  p  £  (Bn  H(t))* ,  H(t)p  and  ext pH(t)p  are  constructed  by  algorithm  III, 
and  uh  1  <  l  <  L,  is  a  complete  system  of  right  H(t)p- coset  representatives  in  H. 

Consider  a  £  KG.  The  component  of  a  in  the  irreducible  G-invariant  subspace  KGt 7 
is  x 

— — -rary,  7  G  extpH(T)p, 

\G{T)p\ 

where  G(t)p  =  A  H (t)p.  Relative  to  the  basis 

uiT'y,  1  <  l  <  L, 

of  KGr-y,  we  can  write 

L 

(XT')  —  aT7(/)u/T7,  ctri(l)  £  K,  7  £  ext pH(t)p. 

(=1 

Theorem  46  For  7  €  ext pH(t)p, 

aT1( l)  =  Y,  oc(xuiw)rl(x~1)\  7(^_1)> 

w£H{t)p  \xeA  / 

where  ti  =  vnruj1 ,  1  <  l  <  L. 

Proof  By  the  results  in  46,  for  every  t  E  G,  we  can  uniquely  write 

t  =  xuiw,  x  £  A,  1  <  l  <  L,  w  £  H(t)p, 

and  we  have 

ir7  =  Ti(x~1)'y(w~1)uiT'). 

Then  v 

0x7  =  £a(()tT7  =  £  E  £  a(xuiw)Tt( x  1)j7(to  *)  J  U/T7, 
f€G  (=1  \w€H{t)p  \x<SA  )  ) 

completing  the  proof. 

Algorithm  V  Expansion  coefficients 

Suppose  a  =  —  Evg//  Eeg.4  oc(xy)xy  is  in  KG. 

A.  For  y  £  H, 

•  Form  ocy  :  A  — »  A  by 

ay(x)  —  a(xy),  x  £  A. 

•  Compute  ay  :  A*  — >  K  by 

&y(T)  =  Y  oiy(x)T(x~1),  r  £  A*. 

xeA 
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linking 

•  Implement  step  A  for  every  y  E  H. 

•  Form  x  :  Y  x  A*  — >  K  by 

x(v,t)  =  3y(r),  y  E  H,t  E  A*. 

B.  For  t  E  A*  and  p  E  {B  fl  #(t))\ 

•  Construct  H(t)p  and  ext pH(t)p  (algorithm  9.3), 

•  Construct  it,,  1  <  l  <  L,  a  complete  system  of  right  H{r)p- coset  representatives 
in  H, 

•  Form  x/  :  H(r)p  — »  K,  1  <  l  <  L,  by 

Xiiw)  =  xiuiTUi'iUtw),  w  E  H(t)p, 

•  Compute  X,  :  ext pH(t)p  — >  AT,  1  <  l  <  L,  by 

Xi(l)=  XiM'riw-1),  7  E  extpH(T)p. 

By  Theorem  46 

ctTJ{l)  =  Xi( 7),  7  €  ext pH{t)p,  1  <  l  <  L. 

C.  Implement  stage  B  for  all  r  E  A*  and  p  E  (B  Pi  H (r))  . 

Stage  C  completes  the  computation  of  the  expansion  coefficient  of  a. 

The  stage  B  computations  of  fc,  1  <  i  <  L,  closely  resemble  Fourier  transform  com¬ 
putations  but  H(t)p  is  not  necessarily  abelian  and  the  computations  are  taken  only  over 
7  G  ext pH(r)p.  We  will  show  how  to  turn  the  main  part  of  these  computations  into  Fourier 
transform  computations. 

Choose  a  complete  system  vjt  1  <  j  <  J,  of  right  {B  C  77(r))-coset  representatives  in 
H{r)p.  Each  w  E  H(t)p  can  be  uniquely  written  as 

w  =  VjU,  1  <  j  <  J,  u  E  B  fl  H(t). 

The  computation  of  x,,  1  <  /  <  L,  can  be  written  as 

X((7)  =  Xi(w)7(w~1) 

w€H(r)p 

=  (vj1)  xiuiTuf1  ,UiVju)p{u~l).  (9) 

j— 1  u^BDH(t) 

Using  this  result  we  can  replace  step  B  by  the  following. 
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B\  For  r  e  A*, 

•  Construct  wk,  1  <  k  <  K ,  a  complete  system  of  right  Bf]H(r )-coset  representa¬ 
tives  in  H , 

•  Form  'tpk  '•  B  fl  H(t )  — >  K,  1  <  k  <  K,by 

ipk(u)  —  xiwki"Wk  1>wk‘u),  uEBC\H(t), 

•  Compute  4>k  ■  (B  fl  H(t))*  — >  K,  1  <  A:  <  K,  by 

i>k{p)  =  52  i>k(u)p(u~1),  p  G  (B  n  H (r))  . 

uGBCiH(t) 


C\  For  pe  (Br\H(r))*, 

•  Construct  ui,  1  <  l  <  L,  a  complete  system  of  right  H(r)p-coset  representatives 
in  H, 

•  Construct  vJt  l<j<J,&  complete  system  of  right  J5nH(r)- coset  representatives 
in  H{r)p. 

linking  For  each  pair  (l,j),  1  <  l  <  L,  1  <  j  <  J,  there  exists  a  unique  k(l,j),  1  < 
k(l,j)  <  K  such  that 

UiVj  =  Wk(ij)Z(i,j),  Z(i,j )  €  Bn  H(t). 

•  Compute  the  products 

$ij{p)  =  p{z(hi))$k{ij)(p),  l<l<L,\<j<J. 


D’.  For  7  G  ext pH(r)p, 


•  Compute  the  matrix  product 


!<(<£,,  l<j<J 


7(«i *) 

livj1) 


OiryO-) 

OCT-f(R) 


By  ( 9)  the  resulting  vector  is 
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5  Code  Development 

Distinct  realizations  of  an  abstract  group  leads  to  a  distinct  unitary  transformations,  and 
distinct  geometric  identification  properties. 

Based  on  the  algorithms  of  section  4,  efficient  codes  implementing  unitary  transforms 
associated  with  several  realizations  of  the  following  abstract  groups  have  been  developed. 

•  (CNxCN)iC2, 

•  (Cjv  x  Cn)  (C*2  x  C2), 

•  (Cjv  X  Cjv)  "3-  {C2  X  C2  X  C2), 

•  (Cjv  X  Cjv)  ( C*2  X  C2  X  C2  X  C2), 

.  (cNxcN)ic3, 

•  (Cjv  X  Cjv)  ^  C4, 

•  (Cjv  X  Cjv)  C6, 

•  (Cjv  x  Cjv)  ^  ((C2  x  C2)  ^  C2), 

•  (Cjv  x  Cjv)  ^  (C4  C2). 

•  (Cjv  <):  C2)  x  (Cm  ^  C2), 

These  codes  lead  to  efficient  implementation  of  the  nonabelian  convolutions. 

Fixing  on  a  particular  representation  of  CN  x  CN,  varying  realizations  of  G  =  ( CN  x 
Cjv)  4  H  is  achieved  by  constructing  2x2  integer  matrices  of  finite  order,  modulo  N.  For 
example,  a  realization  of  the  group  (Cjv  x  CN )  i  (C2  x  C2)  is  achieved  by  constructing 
2x2  integer  matrices  of  order  2  modulo  N,  that  commute.  Viewing  the  elements  of  CN  x 
Cjv  as  column  vectors  of  length  2,  the  group  multiplication  in  G  is  identified  with  matrix 
multiplication. 

Geometric  properties  of  such  realizations  depend  on  the  following. 

•  The  integer  N. 

•  The  integer  matrices  of  finite  order,  modulo  N. 

•  Interaction  between  the  integer  matrices. 
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5.1  Examples 

Denote  the  elements  of  CN  x  CN  by  ambn,  0  <  m,  n  <  N  -  1,  and  characters  of  CN  x  CN 
by  Qkfii,  0  <  k,  l  <  N  —  1, 

TV— i  N  l  Jctyi  *4-  In 

=  E  E 

n=0  m—0 

Denoting  the  generators  of  C2  x  C2  by  s  and  t,  complete  system  of  orthogonal  idempotents 
of  C2  x  C'2  are 

^(1  +  s  +  t  +  st ),  ^(1  —  s  +  t—  st ),  -(1  +  s  —  t  —  st ),  -(1  —  s  —  t  +  st). 

Primitive  idempotents  of  C  ( (CVy  x  Ctv)  {C2  x  C2))  are  of  dimensions  1,  2  or  4.  The  num¬ 
ber,  dimension  and  characteristics  of  primitive  idempotents  depend  on  specific  realizations 

of  C2  xC2. 

Each  realization  is  identified  by  a  group  of  4  commuting  integer  matrices  of  order  2.  Any 
2x2  commuting  matrices  a  and  r  determines  a  realization  of  (Cn  x  Cn)  <$  ( C2  x  C2)  as 

follows.  Set  , 

/v  Ml  «2 
5  *4  O  =  _  , 


t  G4  T  = 


S3  S4  J 

tl  <2 

is  U  ’ 


and  for  0  <  m,  n  <  AT,  let 


a  ,  modiV, 

n 


=r  m  ,  modAT, 
nt  J  [  n  J 

The  realization  of  (CN  x  CV)  <£  (C2  x  C2)  with  a  and  r  is  given  by 

(ambn)a  — »  am’bn’,  {ambn)T  — >  om‘nt. 

The  affine  realizations  are  based  on  the  homogeneous  coordinate  system.  For  a  2  x  2 
matrix  a  of  order  2,  associate  the  matrices 

r  •  n  i 

S  i  52  'I'W 

ffij  =  s3  S4  jj  ,  o  <  i,  j  <  1. 

0  0  1 

i7  ,  is  of  order  2  and  acts  as  an  affine  motion  on  (ro,n)  G  Z/N  x  Z/N  through  the  matrix 
multiplication 

T  m 

cnj  n  ,  0<t,j<l. 
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5.2  Code  development  for  successive  processing 

In  general,  the  results  of  successively  processing  an  image  by  projection  operators  associated 
with  two  distinct  realizations  of  an  abstract  group  are  distinct.  Moreover,  the  order  in  which 
the  projections  are  applied  results  in  significant  differences  in  the  resulting  image. 

One  framework  for  predictable,  successive  processing  is  developed  by  constructing  two 
realizations  A  and  A  <$C  such  that  B  and  C  commute.  In  the  case  of  G  =  ( CN  x 
Cn)  (C2  x  C2),  if  the  two  realizations  of  G  are  obtained  from  4  commuting  matrices 
representing  the  four  copies  of  C2’  s,  then  the  projections  are  not  affected  by  the  order 
of  processing,  and  the  discriminating  characteristics  are  predictable.  The  discriminating 
properties  of  the  resulting  projections  are  equivalent  to  processing  by  the  projection  operators 
associated  with  the  much  larger  group  {CN  x  CN)  4  {@2  x  C2  x  C2  x  C2). 

Example:  Two  commuting  realizations 

The  following  4  elements  commute  and  generate  a  16-group. 


-1  o  1 

[oil 

N  N  I  ' 

0  2 

0  f-1" 

r-H 

i 

o 

1 

> 

i - 

1— ‘ 

o 

1 _ 

? 

to|S: 

1  ' 
l—1 

t 

to|^ 

> 

[f-1  0  J 

All  the  distinct  centralizers  are  realized  by  successively  processing  with  A  =  {CN  x 
CN)  3  (a  x  b)  and  B  —  (Cn  x  Cn)  <£  (c  x  d).  Projections  are  distinguished  by  the  distinct 
subgroups  of  ax  b  and  cx  d.  In  terms  of  imaging  characteristics,  these  projections  distinguish 
the  same  projections  as  the  subgroups  of  the  group  generated  by  all  four  elements.  How¬ 
ever,  the  order  of  the  distinguishing  subgroups  determine  the  dimensionality  of  the  invariant 
subspaces  and  affects  the  detectability  as  opposed  to  discrimination. 

We  have  investigated  and  identified  the  relevant  image  processing  properties  of  many 
realizations  based  on  the  following  MATLAB  codes. 

•  codes  generating  exhaustive  set  of  realizations, 

•  codes  constructing  complete  set  of  primitive  idempotents, 

•  codes  constructing  unitary  transformations  from  complete  set  of  primitive  idempotents, 

•  codes  implementing  varying  filters  in  the  spectral  domain. 
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6  Implementation  results 

All  figures  in  this  report  are  log-scaled,  gray-scale  intensity  plots. 

6.1  Classifying  geometric  properties 

6.1.1  Realizations  of  (Cjv  x  CV)  ^  ( Ci  x  C2) 

Filters  identifying  lines  at  varying  orientations 

Large  collection  of  filters  that  distinguish  lines  at  varying  orientations  have  been  imple¬ 
mented.  For  digital  implementation,  orientations  are  described  by  varying  slopes.  For  the 
selected  realizations  we  chose,  it  is  possible  to  distinguish  between  lines  which  differ  in  slopes 
of  jj.  For  a  large  N,  this  translates  to  small  angular  differences. 
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Filters  composed  of  2-dimensional  idempotents  discriminate  lines  of  slopes  1  and  1. 
An 
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Figure  6.1.1 


Figure  6.1.2 


Figure  6.1.3 


Figure  6.1.4 


Figure  6.1.5 


Figure  6.1.6 


Figure  6.1.1.  Random  intensity  between  0-10  along  lines  of  slope  -1. 

Figure  6.1.2.  Several  lines  of  maximum  intensity  5. 

Figure  6.1.3.  Sum  of  intensities  in  Figures  1.1  and  1.2  which  is  the  input  to  the 
nonabelian  filters. 

Figures  6.1.4  and  6.1.5.  Results  of  filtering  by  2  of  the  2-dimensional  idempotents. 
Figure  6.1.6.  Result  of  filtering  by  the  4-dimensional  idempotents. 
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Performance  of  the  2-dimensional  idempotents  of  C  (( CN  xCN)^(C2x  C2))  in  noise. 


Figure  6.2.1  Figure  6.2.2  Figure  6.2.3 


Figure  6.2.1.  Random  noise  of  variance  15  and  zero  mean. 

Figure  6.2.2.  Several  lines  of  maximum  intensity  5. 

Figure  6.2.3.  Sum  of  intensities  in  Figures  2.1  and  2.2  which  is  the  input  to  the 
nonabelian  filters. 

Figures  6.2.4  and  6.2.5.  Results  of  filtering  by  2  of  the  2-dimensional  idempotents. 
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Filters  composed  of  2-dimensional  idempotents  discriminate  lines  of  slopes  -2  and  0. 


Figure  6.3.4 


Figure  6.3.5 


Figure  6.3.6 


Figure  6.3.1.  Random  intensity  between  0-10  along  lines  of  slope  -1. 

Figure  6.3.2.  Several  lines  of  maximum  intensity  5. 

Figure  6.3.3.  Sum  of  intensities  in  Figures  6.3.1  and  6.3.2  which  is  the  input  to  the 
nonabelian  filters. 

Figures  6.3.4  and  6.3.5.  Results  of  filtering  by  2  of  the  2-dimensional  idempotents. 
Figure  6.3.6.  Result  of  filtering  by  the  4-dimensional  idempotents. 
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Filters  identifying  varying  periodic  structures 
Checkerboard  pattern  is  projected  onto  1-dimensional  idempotents. 


A  = 


0 

-1 


0  1 
1  0 


Figure  6.4.1 


Figure  6.4.1.  32  x  32  alternating  intensities  of  0  and  1. 

Figure  6.4.2.  The  result  of  filtering  by  1-dimensional  idempotents. 


Figure  6.4.3  Figure  6.4.4  Figure  6.4.5 


Figure  6.4.3.  Random  noise  of  variance  1  and  norm  16.19  is  added  to  the  checkerboard 
pattern  in  Figure  6.4.1,  and  used  as  the  input  to  filters  generated  by  the  group  A.  The 
norm  of  the  composite  image  is  32.12. 

Figure  6.4.4.  The  result  of  filtering  by  1-dimensional  idempotents. 

Figure  6.4.5.  The  result  of  filtering  by  4-dimensional  idempotents. 
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Checkerboard  pattern  with  detail  is  projected  into  the  1-dimensional  idempotents  and 
first  of  the  three  2-dimensional  idempotents. 


Figure  6.5.1.  Input  image  of  32  x  32  checkerboard  pattern  is  made  up  of  2  x  2  squares 
of  intensities  2,  3  and  4. 

Figure  6.5.2.  Result  of  filtering  by  1-dimensional  idempotents. 

Figure  6.5.3.  Result  of  filtering  by  the  first  of  the  three  2-dimensional  idempotents. 

Figure  6.5.4 


Figure  6.5.4.  Sum  of  the  intensities  of  images  in  Figures  6.5.2  and  6.5.3. 
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Performance  of  these  filters  are  tested  in  noise. 
Figure  6.5.5  Figure  6.5.6 


Figure  6.5.7 


Figure  6.5.5.  Random  noise  of  values  0-4  is  added  to  the  checkerboard  pattern  of 
Figure  6.5.1  to  form  the  input  image. 

Figure  6.5.6.  The  sum  of  the  results  of  filtering  by  1-dimensional  idempotents  and  the 
first  of  the  three  2-dimensional  idempotents. 

Figure  6.5.7.  The  result  of  filtering  by  the  4-dimensional  idempotents. 
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Checkerboard  pattern  with  detail  is  projected  into  the  1-dimensional  idempotents  and 
second  of  the  three  2-dimensional  idempotents. 


Figure  6.6.1.  Input  image  of  32  x  32  is  made  up  of  4  x  4  squares  of  intensities  of  integer 
values  2-8. 

Figure  6.6.2.  Result  of  filtering  by  1-dimensional  idempotents. 

Figure  6.6.3.  Result  of  filtering  by  the  second  of  the  three  2-dimensional  idempotents. 


Figure  6.6.4  Figure  6.6.5  Figure  6.6.6 


Figure  6.6.4.  Input  image  of  32  x  32  is  generated  by  adding  random  noise  of  intensities 
0-8  to  the  image  in  Figure  6.6.1. 

Figure  6.6.5.  The  sum  of  results  of  filtering  by  1-dimensional  idempotents  and  the 
second  of  the  2-dimensional  idempotents. 

Figure  6.6.6.  Result  of  filtering  by  the  4-dimensional  idempotents. 
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Periodic  structure  is  projected  into  the  second  of  the  three  2-dimensional  idempotents. 
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Figure  6.7.1 


Figure  6.7.1.  Image  of  32  x  32  of  intensity  0-4. 


Figure  6.7.2 


Figure  6.7.3 


Figure  6.7.4 


Figure  6.7.2.  Input  image  of  32  x  32  is  generated  by  adding  random  intensities  of  0-4 
to  the  image  in  Figure  7.1. 

Figure  6.7.3.  The  result  of  filtering  by  the  second  of  the  2-dimensional  idempotents. 
Figure  6.7.4.  Result  of  filtering  by  the  4-dimensional  idempotents. 
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Periodic  pattern  is  projected  into  the  second  of  the  three  2-dimensional  idempotents. 


Figure  6.8.1 


Figure  6.8.1.  Image  of  32  x  32  of  intensity  0-2. 


Figure  6.8.2  Figure  6.8.3  Figure  6.8.4 


Figure  6.8.2.  Input  image  of  32  x  32  is  generated  by  adding  random  noise  of  intensities 
0-4  to  the  image  in  Figure  7.5. 

Figure  6.8.3.  The  result  of  filtering  by  the  second  of  the  2-dimensional  idempotents. 
Figure  6.8.4.  Result  of  filtering  by  the  4-dimensional  idempotents. 

6.1.2  Realizations  of  {Cm  x  Cm)  ^ 

Numerical  experiments  included  in  this  section  illustrate  the  image  processing  characteristics 
of  the  varying  realizations  of  {Cm  x  Cm)  C4.  Robustness  of  these  properties  are  tested  in 
random  white  noise. 
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Figure  6.9.1  Figure  6.9.2 


Figure  6.9.3  Figure  6.9.4 


White  noise  of  zero  mean  and  standard  deviation  of  50%  of  maximum  intensity  is 
added  to  the  image  in  Figure  9.1  to  generate  the  input  image  in  Figure  6.9.2. 

Figure  6.9.3.  The  result  of  projections  by  lower  dimensional  idempotents. 

Figure  6.9.4.  The  residual. 
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Figure  6.10.1  Figure  6.10.2 


Figure  6.10.3  Figure  6.10.4 


White  noise  of  zero  mean  and  standard  deviation  of  maximum  intensity  is  added  to 
the  image  in  Figure  6.10.1  to  generate  the  input  image  in  Figure  6.10.2. 

Figure  6.10.3.  The  result  of  projections  by  lower  dimensional  idempotents. 

Figure  6.10.4.  The  residual. 
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Figure  6.11.1  Figure  6.11.2 


Figure  6.11.3  Figure  6.11.4 


White  noise  of  zero  mean  and  standard  deviation  of  20%  of  the  maximum  intensity  is 
added  to  the  image  in  Figure  11.1  to  generate  the  input  image  in  Figure  6.11.2. 

Figure  6.11.3.  The  result  of  projections  by  lower  dimensional  idempotents. 

Figure  6.11.4.  The  residual. 
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Figure  6.12.4 


White  noise  of  zero  mean  and  standard  deviation  of  50%  of  the  maximum  intensity  is 
added  to  the  image  in  Figure  6.12.1  to  generate  the  input  image  in  Figure  6.12.2. 

Figure  6.12.3.  The  projection  by  lower  dimensional  idempotents. 

White  noise  of  zero  mean  and  standard  deviation  of  the  maximum  intensity  is  added 
to  the  image  in  Figure  6.12.4  to  generate  the  input  image  in  Figure  6.12.5. 

Figure  6.12.6.  The  result  of  projection  by  lower  dimensional  idempotents. 

White  noise  of  zero  mean  and  standard  deviation  of  150%  of  the  maximum  intensity 
is  added  to  the  image  in  Figure  6.12.7  to  generate  the  input  image  in  Figure  6.12.8. 

Figure  6.12.9.  The  projection  by  lower  dimensional  idempotents. 
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6.2  Extending  the  library 

To  extend  the  library  of  projections  discriminating  prescribed  geometry,  we  have  imple¬ 
mented  the  convolution  and  projection  algorithms  associated  with  groups  that  operate  on 
spaces  that  are  unitarily  equivalent  to  the  image  space.  These  include  in  particular,  the 
spaces  of  one  and  two-dimensional  frequencies. 

Figure  6.13.1  Figure  6.13.2  Figure  6.13.3 


White  noise  of  zero  mean  and  standard  deviation  of  100%  of  maximum  intensity  is 
added  to  the  image  of  size  64  x  64  in  Figure  6.13.1  to  generate  the  input  image  in 
Figure  6.13.2. 

Figure  6.13.3  is  the  result  of  nonabelian  group  projection  discriminating  periodicity  in 
frequency. 


Figure  6.14.1 


Figure  6.14.2 


Figure  6.14.3 


White  noise  of  zero  mean  and  standard  deviation  of  300%  of  maximum  intensity  is 
added  to  the  image  of  size  64  x  64  in  Figure  6.14.1  to  generate  the  input  image  in 
Figure  6.14.2. 

Figure  6.14.3  is  the  result  of  nonabelian  group  projection  discriminating  periodicity  in 
frequency. 
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Figure  6.15.1  Figure  6.15.2  Figure  6.15.3 


White  noise  of  zero  mean  and  standard  deviation  of  500%  of  maximum  intensity  is 
added  to  the  image  of  size  64  x  64  in  Figure  6.15.1  to  generate  the  input  image  in 
Figure  6.15.2. 

Figure  6.15.3  is  the  result  of  nonabelian  group  projection  discriminating  periodicity  in 
frequency. 

Figure  6.16.1  Figure  6.16.2  Figure  6.16.3 


White  noise  of  zero  mean  and  standard  deviation  of  100%  of  maximum  intensity  is 
added  to  the  image  of  size  128  x  128  in  Figure  6.16.1  to  generate  the  input  image  in 
Figure  6.16.2. 

Figure  6.16.3  is  the  result  of  nonabelian  group  projection  discriminating  periodicity  in 
frequency. 
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Figure  6.17.1 


Figure  6.17.2 


Figure  6.17.3 
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Figure  6.18.1  Figure  6.18.2  Figure  6.18.3 


White  noise  of  zero  mean  and  standard  deviation  of  500%  of  maximum  intensity  is 
added  to  the  image  of  size  128  x  128  in  Figure  6.18.1  to  generate  the  input  image  in 
Figure  6.18.2. 

Figure  6.18.3  is  the  result  of  nonabelian  group  projection  discriminating  periodicity  in 
frequency. 

6.3  Software  design  tool 

The  library  of  codes  implementing  nonabelian  group  convolutions  and  projections  is  ex¬ 
tensive.  For  image  data  consisting  of  intensities,  the  filters  implemented  in  the  projection 
domain  respects  this  information.  The  result  of  such  filtering  is  an  image  data  of  similar  char¬ 
acteristics,  allowing  successive  and  iterative  processing.  By  combining  several  projections 
and  filters,  extensive  and  non-traditional  image  segmentation  is  possible. 
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6.3.1  Segmentation  by  spatial  periodicity 


Figure  6.19.1 
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Figure  6.19.2 


Figure  6.19.3 


Figure  6.19.1.  Image  of  32  x  32  of  alternating  intensity  of  0  and  2. 

Figure  6.19.2.  Image  of  32  x  32  checkerboard  pattern  made  up  of  2  x  2  squares  of 
intensities  2-5. 

Figure  6.19.3.  The  sum  of  the  intensities  in  Figures  6.19.1  and  6.19.2  is  used  as  input 
for  analysis. 


Figure  6.19.4.  The  result  of  projecting  the  image  in  Figure  6.19.1. 
Figure  6.19.5.  The  sum  of  the  results  of  projecting  the  residual. 
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Figure  6.19.7 


Figure  6.19.7.  Randomized  intensity  of  0-10  added  to  the  image  in  Figure  6.19.3 
Figure  6.19.8.  The  result  of  projecting  the  image  in  Figure  6.19.7. 

Figure  6.19.9.  The  sum  of  the  results  of  projecting  the  residual. 

Figure  6.19.10.  The  residual. 
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Figure  6.20.4 


Figure  6.20.5 


Figure  6.20.1.  Simulated  checkerboard  pattern  of  size  160  x  160. 

Figure  6.20.2.  Downloaded  web  image  of  size  160  x  160. 

Figure  6.20.3.  Sum  of  images  in  Figures  6.20.1  and  6.20.2.  is  used  as  an  input  to 
processing. 

Figures  6.20.4,  6.20.5.  Results  of  projection. 
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Figure  6.21.1 


Downloaded  web  image  of  six  256  x  256. 


Figure  6.21.2 


Figure  6.21.3 


Figure  6.21.4 


Figure  6.21.2.  White  noise  is  added  to  the  image  in  Figure  6.21.1  to  generate  the  input 
image  to  processing. 

Figures  6.21.3  and  6.21.4.  Results  of  projections  discriminating  lines  of  distinct  orien¬ 
tations  and  sizes. 
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Figure  6.21.5 


Figure  6.21.6 


Figure  6.21.7 


Figure  6.21.5.  White  noise  is  added  to  the  image  in  Figure  6.21.1  to  generate  the  input 
image  to  processing. 

Figures  6.21.6  and  6.21.7.  Results  of  projections  discriminating  lines  of  distinct  orien¬ 
tations  and  sizes. 


Figure  6.21.8 


Figure  6.21.9 


Figure  6.21.10 


Figure  6.21.6.  White  noise  is  added  to  the  image  in  Figure  6.21.1  to  generate  the  input 
image  to  processing. 

Figures  6.21.9  and  6.21.10.  Results  of  projections  discriminating  lines  of  distinct  ori¬ 
entations  and  sizes. 
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7  Design  Examples 

7.1  Example:  image  segmentation 

SAR  image  off  the  web  is  segmented  in  two  distinct  ways. 


Figure  7.1.1 
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Download  web  image  of  Stanwick. 
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7.1.1  Image  segmentation  by  spatially  varying  projections 

The  projections  are  constructed  from  the  1-dimensional  idempotents  of  (Cn  x  Cn)  C4. 


Figure  7.1.2 
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Figure  7.1.3 
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Figure  7.1.2.  The  result  of  applying  the  spatially  varying  projection  constructed  from 
1-dimensional  idempotents  of  one  of  the  realizations  of  (Cn  x  Cn)  C±. 


Figure  7.1.3.  The  residual. 
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7.1.2  Image  segmentation  by  textural  differences 

The  projections  are  constructed  from  the  2-dimensional  idempotents  of  (CW  x  Cn)  <£  C\. 


Figure  7.1.4 
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Figure  7.1.5 
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Figure  7.1.4.  The  result  of  applying  the  projection  discriminating  textural  (periodic) 
information  constructed  from  2-dimensional  idempotents  of  one  of  the  realizations  of 
(< CN  x  CN)  4  C4. 

Figure  7.1.5.  The  residual. 
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Sidescan  sonar  image  is  segmented  by  textural  differences. 

Figure  7.1.6 


Sidescan  sonar  image  of  ocean  floor. 
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Figure  7.1.7 
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Figure  7.1.8 
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Figure  7.1.7.  The  result  of  applying  the  projection  discriminating  textural  (periodic) 
information  constructed  from  2-dimensional  idempotents  of  one  of  the  realizations  of 
(CNXCN)4CA. 


Figure  7.1.8.  The  residual. 
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7.2  Software  tool  for  automatic  pattern  localization 

Codes  implementing  detection  of  prescribed  patterns  have  been  developed  to  locate  the 
position  of  the  patterns.  Unlike  the  traditional  methods  based  on  matched  filtering,  the 
nonabelian  projection  methods  can  detect  and  locate  patterns  of  abstract  descriptions  that 
are  independent  of  variances  in  scale  and  magnitude  of  the  patterns. 

Encoding  of  prescribed  pattern  into  an  invariant  subspace  is  a  design  process,  based  on 
selection  of  a  nonabelian  group.  Once  this  encoding  is  completed  however,  locating  the 
pattern  is  automatic. 

A  pattern  that  can  be  encoded  into  an  invariant  subspace  is  embedded  into  various 
backgrounds.  The  nonabelian  group  chosen  for  the  example  is  one  realization  of  the  abstract 
group  G  =  (Cjv  x  Cjv)  $  (C2  x  C2), 


C2  x  C2  ~ 


The  operator  which  projects  an  image  into  the  selected  invariant  subspace  partitions  the 
image  within  the  projection  (depending  on  the  location  of  the  encoded  pattern  and  the 
background)  into  quadrants,  by  4  levels  of  magnitudes,  by  the  action  of  C2  x  C2-  We  have 
used  the  relative  levels  of  magnitudes  of  the  quadrants  to  determine  the  amount  of  the 
prescribed  pattern  occupying  the  quadrants.  From  the  knowledge  of  the  relative  size  of 
the  prescribed  pattern,  off-set  location  is  determined  successively.  Current  implementation 
requires  4  iterations.  The  number  of  iterations  depend  largely  on  the  extent  of  the  pattern: 
The  larger  the  pattern,  the  fewer  the  iteration.  At  each  iteration,  the  test  data  size  is  reduced 
by  50%.  The  computational  complexity  of  each  iteration  is  in  the  order  of  2AflogAT,  where 
N  is  the  test  data  size. 

The  current  implementation  is  parameterized  for  N,  a  power  of  2,  which  determines  the 
data  size,  2 N  x  2N  and  the  size  of  the  pattern  is  set  at  y  x  y. 
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7.2.1  Example  scenes 

Location  of  horizontal  bars  in  random  background 


Figure  7.2.4 


Figure  7.2.5 


Figure  7.2.6 


Figure  7.2.1.  Pattern  of  alternating  horizontal  bars  of  size  16  x  16. 

Figure  7.2.2.  Pattern  in  Figure  3.1  is  placed  in  the  scene  of  randomized  magnitude  of 
size  64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.2.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (0, 0),  the  test  image  on  coordinates  [0  :  47, 0  :  47]  is  displayed. 

Figure  7.2.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (0, 16),  the  test  image  on  coordinates  [0  :  31, 16  :  47]  is  displayed  . 

Figure  7.2.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (8, 16),  the  test  image  on  coordinates  [8  :  31, 16  :  39]  is  displayed  . 

Figure  7.2.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (10,20),  the  test  image  on  coordinates  [10  :  29,20  :  39]  is  displayed. 
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Location  of  noisy  horizontal  bars  in  random  background 


Figure  7.3.1  Figure  7.3.2  Figure  7.3.3 


Figure  7.3.4  Figure  7.3.5  Figure  7.3.6 


Figure  7.3.1.  Pattern  of  alternating  noisy  horizontal  bars  of  size  16  x  16. 

Figure  7.3.2.  Pattern  in  Figure  4.1  is  placed  in  the  scene  of  randomized  magnitude  of 
size  64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.3.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (0, 16),  the  test  image  on  coordinates  [0  :  47, 16  :  63]  is  displayed. 

Figure  7.3.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (0,  32),  the  test  image  on  coordinates  [0  :  31,32  :  63]  is  displayed  . 

Figure  7.3.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (8,36),  the  test  image  on  coordinates  [8  :  31,36  :  59]  is  displayed  . 

Figure  7.3.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (12,38),  the  test  image  on  coordinates  [12  :  31,  38  :  57]  is  displayed. 
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Location  of  horizontal  bars  in  random  background 


Figure  7.4.1  Figure  7.4.2  Figure  7.4.3 


Figure  7.4.4  Figure  7.4.5  Figure  7.4.6 


Figure  7.4.1.  Pattern  of  horizontal  bars  of  size  16  x  16. 

Figure  7.4.2.  Pattern  in  Figure  5.1  is  placed  in  the  scene  of  randomized  magnitude  of 
size  64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.4.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (16, 16),  the  test  image  on  coordinates  [16  :  63, 16  :  63]  is  displayed. 

Figure  7.4.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (32,32),  the  test  image  on  coordinates  [32  :  63,32  :  63]  is  displayed  . 

Figure  7.4.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (36, 32),  the  test  image  on  coordinates  [36  :  59, 32  :  55]  is  displayed  . 

Figure  7.4.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (42, 39),  the  test  image  on  coordinates  [42  :  61, 39  :  58]  is  displayed. 
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Location  of  noisy  horizontal  bars  in  random  background 


Figure  7.5.1  Figure  7.5.2  Figure  7.5.3 


Figure  7.5.4  Figure  7.5.5  Figure  7.5.6 


Figure  7.5.1.  Pattern  of  noisy  horizontal  bars  of  size  16  x  16. 

Figure  7.5.2.  Pattern  in  Figure  6.1  is  placed  in  the  scene  of  randomized  magnitude  of 
size  64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.5.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (16, 16),  the  test  image  on  coordinates  [16  :  63, 16  :  63]  is  displayed. 

Figure  7.5.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (16, 24),  the  test  image  on  coordinates  [16  :  47,  24  :  55]  is  displayed  . 

Figure  7.5.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (24, 32),  the  test  image  on  coordinates  [24  :  47, 32  :  55]  is  displayed  . 

Figure  7.5.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (24, 32),  the  test  image  on  coordinates  [24  :  43, 32  :  51]  is  displayed. 


Figure  7.6.1.  Pattern  of  noisy  horizontal  bars  of  size  16  x  16. 

Figure  7.6.2.  Pattern  in  Figure  7.1  is  placed  in  the  scene  of  randomized  magnitude  of 
size  64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.6.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (8, 0),  the  test  image  on  coordinates  [8  :  55, 0  :  47]  is  displayed. 

Figure  7.6.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (16,0),  the  test  image  on  coordinates  [16  :  47, 0  :  31]  is  displayed  . 

Figure  7.6.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (24,0),  the  test  image  on  coordinates  [24  :  47,  0  :  23]  is  displayed  . 

Figure  7.6.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (24,4),  the  test  image  on  coordinates  [24  :  43, 4  :  23]  is  displayed. 
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Location  of  horizontal  bars  in  pseudo-random  background 


Figure  7.7.1  Figure  7.7.2  Figure  7.7.3 


Figure  7.7.1.  Pattern  of  horizontal  bars  of  size  16  x  16. 

Figure  7.7.2.  Pattern  in  Figure  7.71  is  placed  in  the  somewhat  structured  scene  of  size 
64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.7.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (0, 0),  the  test  image  on  coordinates  [0  :  47, 0  :  47]  is  displayed. 

Figure  7.7.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (16, 16),  the  test  image  on  coordinates  [16  :  47, 16  :  47]  is  displayed. 

Figure  7.7.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (16, 16),  the  test  image  on  coordinates  [16  :  39, 16  :  39]  is  displayed. 

Figure  7.7.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (20, 18),  the  test  image  on  coordinates  [20  :  39, 18  :  37]  is  displayed. 
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Location  of  horizontal  bars  in  structured  background 


Figure  7.8.1 


Figure  7.8.2 


Figure  7.8.3 
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Figure  7.8.4 


Figure  7.8.5 


Figure  7.8.6 


Figure  7.8.1.  Pattern  of  horizontal  bars  of  size  16  x  16. 

Figure  7.8.2.  Pattern  in  Figure  7.8.1  is  placed  in  the  somewhat  structured  scene  of  size 
64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.8.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (16,0),  the  test  image  on  coordinates  [16  :  63, 0  :  47]  is  displayed. 

Figure  7.8.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (32, 0),  the  test  image  on  coordinates  [32  :  63, 0  :  31]  is  displayed  . 

Figure  7.8.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (40, 4),  the  test  image  on  coordinates  [40  :  63, 4  :  27]  is  displayed  . 

Figure  7.8.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (42, 6),  the  test  image  on  coordinates  [42  :  63, 6  :  25]  is  displayed. 
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Location  of  horizontal  bars  in  pseudo-structured  background 


Figure  7.9.1 


Figure  7.9.2 


Figure  7.9.3 


Figure  7.9.4 


Figure  7.9.5 


Figure  7.9.6 


Figure  7.9.1.  Pattern  of  horizontal  bars  of  size  16  x  16. 

Figure  7.9.2.  Pattern  in  Figure  7.9.1  is  placed  in  the  somewhat  structured  scene  of  size 
64  x  64.  The  placement  is  chosen  randomly. 

Figure  7.9.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (0, 0),  the  test  image  on  coordinates  [0  :  47, 0  :  47]  is  displayed. 

Figure  7.9.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (0,0),  the  test  image  on  coordinates  [0  :  31,0  :  31]  is  displayed  . 

Figure  7.9.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (8, 0),  the  test  image  on  coordinates  [8  :  31, 0  :  23]  is  displayed. 

Figure  7.9.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (8,0),  the  test  image  on  coordinates  [8  :  27, 0  :  19]  is  displayed. 
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Location  of  horizontal  bars  in  pseudo-structured  background 


Figure  7.10.1 


Figure  7.10.2 


Figure  7.10.3 


Figure  7.10.4 


Figure  7.10.5 


Figure  7.10.6 


Figure  7.10.1.  Pattern  of  horizontal  bars  of  size  32  x  32. 

Figure  7.10.2.  Pattern  in  Figure  7.10.1  is  placed  in  the  somewhat  structured  scene  of 
size  128  x  128.  The  placement  is  chosen  randomly. 

Figure  7.10.3.  Based  on  the  result  of  the  first  iteration  of  the  projection  of  off-set 
coordinate  of  (16,32),  the  test  image  on  coordinates  [16  :  111,32  :  127]  is  displayed. 

Figure  7.10.4.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (32,64),  the  test  image  on  coordinates  [32  :  95,  64  :  127]  is  displayed  . 

Figure  7.10.5.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (32,80),  the  test  image  on  coordinates  [32  :  79,80  :  127]  is  displayed  . 

Figure  7.10.6.  Based  on  the  result  of  the  second  iteration  of  the  projection  of  off-set 
coordinate  of  (41,81),  the  test  image  on  coordinates  [41  :  80,81  :  120]  is  displayed. 
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7.3  Software  for  automatic  location  of  targets 

7.3.1  The  problem 

Dr.  Albanese  (210-536-5710)  at  Brooks  AFB  has  posed  the  problem  of  detecting  and  locating 
rectangles  in  severe  noise:  the  level  of  noise  is  such  that  the  rectangles  are  hardly  recognizable. 

Several  of  the  previously  implemented  nonabelian  group  filters  project  line  segments  in 
somewhat  noisy,  cluttered  background.  Digitally  recorded  rectangles  are  viewed  as  several 
line  segments  of  the  same  lengths  occupying  some  consecutive  pixels  in  the  direction  perpen¬ 
dicular  to  the  line  segments.  As  the  noise  level  increases,  and  the  size  of  the  object  reduces, 
the  projections  become  less  robust. 

7.3.2  A  solution 
Imaging  model 

In  general,  an  imaging  model  for  an  image  of  size  N  x  N,  N  =  2K ,  is  based  on 

•  A  semi-direct  product  group  A  $  B  of  order  N2  where  A  and  B  are  abelian  groups, 

•  An  indexing  of  the  image  by  the  group. 

For  reasons  explained  more  fully  in  the  next  section,  we  choose  A  and  B  so  that 

•  A  —  Z/2m  x  Z/2m,  M  =  K  -  2. 

•  B  is  a  subgroup  of  GL{ 2, 2M)  isomorphic  to 

C2  x  C2  x  C2  x  C2, 

where  C2  is  the  cyclic  group  of  order  2. 

We  have  previously  (Report  3,  January  2001)  investigated  various  actions  of  GL( 2, 2M)  on 
A.  To  construct  subgroups  B  of  this  form,  we  need  to  find  four  commuting  elements  of  order 
2  in  GL(2, 2M). 

By  an  action  of  B  on  A,  we  always  mean  a  semi-direct  product  action.  For  a  fixed  action 
of  B  on  A,  the  indexing  of  the  image  by  A  B  is  given  as  in  the  table. 


b0A 

hA 

M 

M 

b5A 

h  A 

bjA 

b8A 

M 

bw  A 

bnA 

h2A 

bis  A 

buA 

bi^A 

Table  1.  Ordering  of  A  B 

The  indexing  of  B  =  {1  =  60,  &i, . . . ,  6i5}  will  have  no  effect  on  the  resulting  processing. 

The  design  problem  is  to  construct  an  action  of  B  on  A  such  that  relative  to  the  specified 
indexing  some  translation  invariant  Subspace  in  the  group  algebra  C(A  ■=)  B)  will  contain 
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lines  in  the  image  at  a  prescribed  orientation.  From  the  general  theory,  the  solution  to  this 
problem  is  given  by  studying  the  induced  action  of  B  on  the  set  of  characters  A*  viewed  as 
a  subset  of  C(A  B).  First  index  A  by  (kuk2),  0  <  ku  h  <  2M  and  consider  A  as  ordered 
lexicographically.  The  set  of  characters  A*  is  indexed  by  {h,h)>  0  <  l\,  l2  <  2M,  where  the 
character  a /,  /,  is  defined  by 

n  .  1 1  k]  -plo  ^2 

ahh(kuk2)  =  e2™^^ 


A*  is  ordered  lexicographically  as  well. 

Suppose  the  design  problem  is  to  choose  an  action  of  B  on  A  such  that  a  line  of  slope 
— c  in  the  image  is  contained  in  some  translation  invariant  subspace  associated  to  A  <$  B  in 
C(A  <£  B).  To  solve  this  problem  we  must  choose  an  action  of  B  on  A  such  that  the  induced 
action  of  B  on  A*  satisfies  the  property  that  there  exists  a  b  G  B  such  that  in  C(A  <1  B), 


0  <  Zx  <  2m. 


This  is  equivalent  to  the  assertion  that  the  centralizer  subgroup  in  C(A  <£  B)  of  the  line  of 
characters,  aduil,  0  <  h  <  2M,  is  nontrivial.  For  the  prescribed  orientation  of  the  lines  in 
the  image  of  the  project  described  above,  the  developed  software  implements  actions  of  B 
on  A  solving  the  design  problem. 

This  software  has  another  feature  which  we  have  found  is  especially  important.  There 
may  be  many  actions  of  B  on  A  which  solve  the  design  problem.  For  producing  the  most 
robust  (in  noise  or  clutter)  filtering  of  images  for  detecting  of  lines  at  prescribed  orientations, 
we  have  found  that  some  are  better  than  others  based  on  the  following  criterion. 

Fix  an  action  6  of  B  on  A  and  form  the  resulting  semi-direct  product  A  eB.  For  a 
line  of  characters  and  a  subgroup  C  of  B,  let  n(C,9)  be  number  of  characters  in  the  line 
centralized  by  C  in  C (A  0B).  As  C  runs  over  the  subgroups  of  B,  the  numbers  n(C,9) 
may  vary.  Denote  the  largest  by  n(9).  As  6  runs  over  the  actions  of  B  on  A ,  the  numbers 
n{9)  may  vary.  We  have  found  that  actions  6  of  B  on  A  producing  the  largest  n(9)  will  yield 
the  most  robust  (in  noise  or  clutter)  detection  of  lines  of  slope  —  c. 


Choice  of  the  order  of  A 

In  general  the  relative  orders  of  A  and  B  in  the  indexing  of  an  image  by  A  $  B  affects  the 
localization  properties  of  the  resulting  nonabelian  group  filtering  operations.  There  is  an 
ambiguity  of  location  determined  by  the  size  of  A,  the  smaller  the  order  of  A, the  smaller  the 
ambiguity.  In  the  ideal  case,  the  size  of  the  object  should  be  the  same  as  the  order  of  A. 
This  has  the  following  consequence  for  the  specific  design  problem  considered  in  the  project 
describe  above. 

The  line  defined  by  a  pair  (-I2J1)  in  A  is  an  algebraic  line  of  length  2M.  An  image 
containing  a  geometric  line  segment  of  the  same  slope  will  have  a  large  coefficient  in  its 
orthogonal  idempotent  expansion.  Although  this  property  completely  localizes  the  geometric 
line  segment  in  the  direction  (h,h),  there  are  ambiguities  in  determining  the  length  and 
position  of  the  line  segment  upto  the  order  of  A.  A  should  be  chosen  so  that  the  lines  in  A 
have  the  same  order  as  the  geometric  line  segment:  ambiguity  in  a  few  pixels  can  be  resolved 


AF98T009  Phase  II:  Design  Example 


105 


by  other  criteria.  In  noise-free  images,  the  relative  magnitudes  of  the  expansion  coefficients 
exactly  determine  the  length  of  a  line  segment. 


Software  tool  for  detection/localization 

The  software  tool  incorporates  projection  operators  from  three  different  realizations  of 

C2  x  C2  x  C2  x  C2, 


to  mitigate  clutter  and  noise.  One  such  realization  is  given  as  the  group  generated  by 


o 

rH 

1 

l _ 

o 

i--4 

1 

1  _ 

O 

to 

S! 

_ i 

ij 

5 

T— 1 

i 

o 

_ 1 

> 

2m  x 

> 

2m  -  1  0 
2m  2m  -  1 


Each  of  the  three  realizations  of  B  contains  the  element 


b  = 


-1 

-1 


0 

1 


Observe  that  b  £  B(a o^),  0  <  l\  <  2M . 

The  software  tool  processes  an  image  of  size  128  x  128  by  projecting  the  image  into 
orthogonal  invariant  subspaces  of  C(A  <£  B),  using  the  three  different  realizations  of  B.  The 
location  of  the  maximal  coefficients  in  the  spaces  spanned  by 

raduh,  0  <  h  <  2m,  t  £  B(adl:h)* 

are  compared.  If  the  locations  of  the  maximal  coefficients  from  the  three  distinct  projections 
are  the  same,  and  the  y-coordinate  is  produced  and  the  confidence  level  is  set  to  3  (out  of 
0,  1,  2,  and  3).  This  is  due  to  the  fact  that  the  three  distinct  projections  behave  differently 
in  the  presence  of  randomized  white  noise.  If  such  locations  differ  by  less  than  the  width 
of  the  object  of  interest,  the  weighted-average  location  is  produced  as  the  y-coordinate  and 
confidence  level  is  reduced  from  3  by  1  or  2  depending  on  the  values  and  the  locations  of  the 
maximal  coefficients. 

To  localize  the  ^-coordinate  as  well  as  the  length  of  the  object,  parameterized  thresholding 
of  the  projection  coefficients  is  used.  Again,  depending  the  number  of  matches  in  the  three 
different  projections,  ^-coordinate,  the  length  of  the  object  and  updated  confidence  level  is 
produced. 


7.3.3  Example  scenes 

All  the  figures  in  this  report  are  log-scaled,  gray  scale,  intensity  plots.  In  each  set  of  three 
images,  the  image  on  top  left  is  the  noise-free  scene  consisting  of  two  rectangular  boxes  in 
structured  background.  The  rectangular  boxes  vary  in  sizes  up  to  4  pixels  in  each  direction, 
and  placed  at  random  in  the  scene.  The  intensity  of  the  structured  background  is  120%  of 
the  intensity  of  the  boxes.  The  image  on  top  right  is  the  result  of  adding  white  noise  to 
the  scene  on  top  left  and  is  used  as  input  to  processing.  The  software  tool  automatically 
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produces  the  coordinates  of  the  centers  of  the  two  boxes  along  with  a  confidence  level.  The 
first  component  of  the  ordered  pair  indicating  the  coordinate  is  the  vertical  distance  from  the 
top  of  the  image,  while  the  second  component  is  the  horizontal  distance  from  the  left.  The 
image  on  bottom  left  is  the  noise-free  scene  with  the  resulting  coordinates  highlighted.  The 
purpose  of  this  is  to  show  the  positions  of  the  computed  coordinates  relative  to  the  original 
scene. 


20  40  60  80  100  120 


vertical  box,  centered  at  (25,19),  c-level=  3 
horizontal  box,  centered  at  (11,43),  c-level=  3 


85,104, 3;99, 105,3 


ft 
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20  40  60  80  100  120 


vertical  box,  centered  at  (85,104),  c-level=  3 
horizontal  box,  centered  at  (99,105),  c-level=  3 
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20  40  60  80  100  120 


vertical  box,  centered  at  (50,13),  c-level=  3 
horizontal  box,  centered  at  (7,105),  c-level=  3 
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20  40  60  80  100  120 


vertical  box,  centered  at  (83,20),  c-level=  1 
horizontal  box,  centered  at  (99,18),  c-level=  2 
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20  40  60  80  100  120 


vertical  box,  centered  at  (20,40),  c-level=  3 
horizontal  box,  centered  at  (116,105),  c-level=  1 
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20  40  60  80  100  120 


vertical  box,  centered  at  (46,95),  c-level=  3 
horizontal  box,  centered  at  (55,105),  c-level=  1 
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20  40  60  80  100  120 


vertical  box,  centered  at  (105,34),  c-level=  3 
horizontal  box,  centered  at  (106,25),  c-level=  0 
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20  40  60  80  100  120 


vertical  box,  centered  at  (45,77),  c-level=  0 
horizontal  box,  centered  at  (118,64),  c-level=  3 
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20  40  60  80  100  120 


vertical  box,  centered  at  (105,29),  c-level=  3 
horizontal  box,  centered  at  (58,22),  c-level=  3 
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20  40  60  80  100  120 


vertical  box,  centered  at  (51,95),  c-level=  3 
horizontal  box,  centered  at  (26,82),  c-level=  1 
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20  40  60  80  100  120 


thick  box,  centered  at  (24,95),  c-level=  3 
thin  box,  centered  at  (63,83),  c-level=  3 
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20  40  60  80  100  120 


thick  box,  centered  at  (105,87),  c-level=  3 
thin  box,  centered  at  (80,80),  c-level=  3 
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